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ABSTRACT 

We solve the general-relativistic steady-state eigenvalue problem of neutrino-driven protoneutron star winds, 
which immediately follow core-collapse supernova explosions. We provide velocity, density, temperature, and 
composition profiles and explore the systematics and structures generic to such a wind for a variety of protoneu- 
tron star characteristics. Furthermore, we derive the entropy, dynamical timescale, and neutron-to-seed ratio in the 
general relativistic framework essential in assessing this site as a candidate for r-process nucleosynthesis. Gener- 
ally, we find that for a given mass outflow rate (M), the dynamical timescale of the wind is significantly shorter 
than previously thought. 

We argue against the existence or viability of a high entropy (> 300 per ke per baryon), long dynamical 
timescale r-process epoch. In support of this conclusion, we model the protoneutron star cooling phase, cal- 
culate nucleosynthetic yields in our steady-state profiles, and estimate the integrated mass loss. We find that 
transonic winds enter a high entropy phase only with very low M (< 1 X 10~ 9 M Q s" 1 ) and extremely long dynam- 
ical timescale (r p > 0.5 seconds). Our results support the possible existence of an early r-process epoch at modest 
entropy (~ 150) and very short dynamical timescale, consistent in our calculations with a very massive or very 
compact protoneutron star that contracts rapidly after the preceding supernova. We explore possible modifications 
to our models, which might yield significant r-process nucleosynthesis generically. 

Finally, we speculate on the effect of fallback and shocks on both the wind physics and nucleosynthesis. We 
find that a termination or reverse shock in the wind, but exterior to the wind sonic point, may have important nu- 
cleosynthetic consequences. The potential for the r-process in protoneutron star winds remains an open question. 

Subject headings: Winds, neutrinos, supernovae, neutron stars, nucleosynthesis 



1. INTRODUCTION 



The successful two-dimensional Type-II supernova simula- 
tion of Burrows, Hayes, and Fryxell (1995) shows clearly a 
post-explosion neutrino-driven wind, emerging approximately 
half a second after bounce. The convective plumes and fingers 
due to Rayleigh-Taylor instabilities that accompany shock re- 
ignition in the gain region are pushed out and cleared from the 
area closest to the neutron star by the pressure of the neutrino- 
driven wind. The last 50 milliseconds (ms) of the simulation 
show that a nearly spherically symmetric wind has established 
itself as the protoneutron star, newly born, begins its Kelvin- 
Helmholtz cooling phase. 1 Although these simulations em- 
ployed only crude neutrino transport, did not address the issue 
of fallback, and did not study the wind as a function of progen- 
itor mass and structure, they are suggestive of a general phe- 
nomenon that might naturally accompany many core-collapse 
supernovae. 

1 See these 2D simulations at http://www.astrophysics.arizona.edu/movies.html. 



Some multiple of 10 erg will be lost via neutrino radiation 
by the protoneutron star as it cools. A small fraction of that 
energy will be deposited in the surface layers of the nascent 
neutron star, heating and driving material from its surface. Al- 
though the wind is interesting in its own right, hydrodynami- 
cally and as a phenomenon that attends both the supernova and 
the Kelvin-Helmholtz cooling phase, perhaps its most impor- 
tant ramification is the potential production of ^50% of all the 
nuclides above the iron group in rapid(r) neutron-capture nu- 
cleosynthesis. 

In the r-process, rapid interaction of neutrons with 
heavy, neutron-rich, seed nuclei allows a neutron capture- 
disintegration equilibrium to establish itself among the isotopes 
of each element. Beta decays occur on a longer timescale and 
increase the nuclear charge. For sufficiently large neutron-to- 
seed ratio, the 'nuclear flow' proceeds to the heaviest nuclei, 
forming abundance peaks at A ~ 80, 130, and 195 (Burbidge 
et al. 1957; Wallerstein et al. 1997). The neutron-to-seed ra- 
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tio itself is largely set by the dynamical timescale (r p ; see 
eq. 32), entropy (s), and neutron richness in the earlier phase 
of the expansion (Hoffman, Woosley, and Qian 1997; Meyer 
and Brown 1997; Freiburghaus et al. 1999). While the nuclear 
physics is fairly well understood, the astrophysical site has not 
been unambiguously established. Currently, neutron star merg- 
ers (Freiburghaus, Rosswog, and Thielemann 1999, and refer- 
ences therein; Rosswog et al. 1999) and protoneutron star winds 
(Meyer et al. 1992; Woosley and Hoffman 1992; Woosley et 
al. 1994; Qian and Woosley 1996; Hoffman, Woosley, and Qian 
1997; Otsuki et al. 2000) are considered the most viable candi- 
dates. 

In addition to attaining the requisite neutron-to-seed ratios, 
dynamical timescales, and temperatures, the astrophysical site 
must consistently reproduce the observed solar abundance pat- 
tern of r-process elements with A > 135. Recent observations 
of neutron-capture elements in ultra-metal-poor halo stars (Bur- 
ris et al. 2000; McWilliam et al. 1995a,b; Sneden et al. 1996; 
Cowan et al. 1996; Westin et al. 2000; Hill et al. 2001) show re- 
markable agreement with the solar r-process abundance pattern 
in this mass range. This suggests a universal mechanism for 
producing the second and third abundance peaks, which acts 
early in the chemical enrichment history of the galaxy. In this 
paper, beyond addressing the physical nature and systematics 
of the wind, we investigate its potential as a site for r-process 
nucleosynthesis up to and beyond the third abundance peak. 

Duncan, Shapiro, & Wasserman (1986) were the first to ad- 
dress the physics of steady-state neutrino-driven neutron star 
winds. Although interested in the relative importance of the 
neutrino and photon luminosity in determining the wind dy- 
namics, they also identified some of the basic systematics and 
scaling relations generic to the problem. More recent investi- 
gations have focused less on the general physics of the wind 
and more on its potential nucleosynthetic yield. Woosley et 
al. (1994) conducted the first such detailed study. They fol- 
lowed the nucleosynthesis in a protoneutron star wind that 
emerged in a one-dimensional post-supernova environment. 
Approximately 18 seconds after collapse and explosion, their 
model attained entropies of ~400 (throughout, we quote en- 
tropy per ka per baryon), long dynamical timescales, and elec- 
tron fraction (Y e ) in the range 0.36-0.44. However, in their 
model the supernova shock reached only 50,000 km at these 
late times. In turn, this external boundary caused the wind ma- 
terial to move slowly. It remained in the heating regime for an 
extended period, thus raising the entropy above what any simu- 
lation or analytical calculation has since obtained. Although the 
r-process proceeded to the third abundance peak in their calcu- 
lation, nuclei in the mass range near A <~ 90 (particularly, 88 Sr, 
89 Y, and 90 Zr) were overproduced by more than a factor of 100. 

Takahashi, Witti, & Janka (1994) conducted a similar in- 
vestigation, but did not attain the entropies of Woosley et 
al. (1994). In fact, they fell short by a factor of <~ 5. Later, 
Qian & Woosley (1996) showed that for reasonable protoneu- 
tron star characteristics, including post-Newtonian corrections, 
an entropy of 400 is unrealistic. Qian & Woosley (1996) 
also provided many analytical scaling relations that have since 
framed the discussion of neutrino-driven winds. Following this 
work, Hoffman, Woosley, & Qian (1997) conducted nucle- 
osynthetic calculations in the wind models of Qian & Woosley 
(1996). They concluded that the standard wind models of Qian 
& Woosley (1996) did not produce third peak r-process nu- 



cleosynthesis. They also employed an adiabatic cooling pre- 
scription to survey the parameter space relevant to protoneu- 
tron star winds and noted several important systematics in the 
nucleosynthesis. Particularly, they identified a low entropy 
(<~ 120), fast timescale (t p <~ 2 ms), and high electron fraction 
(Y e > 0.48) window where third peak r-process could take place 
(see their Fig. 10). 

Cardall & Fuller (1997) generalized some of the scaling re- 
lations presented in Qian & Woosley (1996) to include general 
relativity and found significant enhancements in the entropy 
and dynamical timescale of the wind in this framework. Re- 
cently, Otsuki et al. (2000) have sought to solve the general- 
relativistic wind equations and to conduct r-process calcula- 
tions in the winds they obtained. They concluded that r-process 
nucleosynthesis can proceed to the third abundance peak at 
A ~ 195 for a protoneutron star with radius 10 km, a gravi- 
tational mass of 2 M Q , and total neutrino luminosity of 10 52 
erg s" 1 . These conditions produce modest entropies (s <~ 140) 
and fast dynamical timescales (r p <~ 2-3 ms). Although a sig- 
nificant parameter in protoneutron star winds, Y e was fixed by 
hand in their calculations. Furthermore, they employed a sim- 
ple equation of state and, in the context of the transonic wind 
problem, an unphysical external boundary condition. 

These uncertainties and ambiguities in the conclusions of 
previous groups suggest that a re-evaluation is in order. Our 
goal in this paper is to solve the full eigenvalue problem of 
the steady-state transonic wind problem in general relativity, 
employing physical boundary conditions. Using this formal- 
ism, we survey the relevant parameter space, identify the major 
systematic trends, and explore some of the particulars of the 
general-relativistic treatment. We then use these steady-state 
solutions to model the whole of the Kelvin-Helmholtz cool- 
ing phase, including radial contraction of the protoneutron star, 
as well as the evolution of the neutrino luminosity and aver- 
age neutrino energy. We estimate the total amount of mate- 
rial ejected during this cooling epoch, and put significant con- 
straints on the range of entropies and dynamical timescales that 
might actually occur in Nature. In addition, for a subset of the 
models generated, we calculate the total nucleosynthetic yield 
as a function of atomic mass. 

In §2, we present the fundamental equations for a time- 
independent neutrino-driven wind (general-relativistic and 
Newtonian), including the equation for the evolution of the 
electron fraction. In addition, we present the integrals of the 
motion and a discussion of the equation of state we employ. In 
§3, we describe our solution to the wind problem using an iter- 
ative relaxation procedure on an adjustable radial grid and the 
necessary boundary conditions. In §4, we present the neutrino 
heating rates used in this study, review the effects of general 
relativity, discuss potential modifications to the energy depo- 
sition rates, and explore (approximately) the effects of trans- 
port. In §5, we present our results for the wind problem itself. 
The wind structures and general characteristics as a function of 
neutron star mass, radius, and neutrino spectral character are 
explored. We cover the entire relevant parameter space and in- 
clude some modifications to the power laws presented in Qian 
& Woosley (1996). In §6, we use our steady-state wind re- 
sults to construct a sequence of such models that represent the 
time evolution of the wind during the protoneutron star cooling 
phase. We estimate the total mass ejected for a given evolution- 
ary trajectory and put useful constraints on possible epochs of 
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r-process nucleosynthesis. In §7, we present nucleosynthetic 
results from a subset of our wind trajectories. In §8, we discuss 
reasonable modifications to our wind models that might yield 
a successful r-process, including changes to the energy depo- 
sition function. Finally, in §9 we review our results, summa- 
rizing the constraints our calculations impose on the viability 
of the protoneutron star wind as the astrophysical site of the r- 
process. Furthermore, we speculate on the effects of progenitor 
structures and fallback, as well as hydrodynamical and trans- 
port considerations left to be addressed in future work. 

2. THE STEADY-STATE WIND EQUATIONS 

The time-independent hydrodynamical equations for flow in 
a Schwarzschild spacetime can be written in the form (Nobili, 
Turolla, and Zampieri 1991; Flammang 1982) 



1 d(yy) | Idp | 2_ Q 
vy dr p dr r 



1 dy 



1 dP 



and 



y dr e+P dr 
de e + P dp 



= 0, 



dr 



dr vy 



(1) 



(2) 



(3) 



where u r (= vy) is the radial component of the fluid four- velocity, 
v is the velocity of the matter measured by a stationary observer, 



y = 



l-2GM/rc 2 



1/2 



(4) 



l-v 2 /c 2 

e (= pc 2 + pe) is the total mass-energy density, p is the rest-mass 
density, P is the pressure, e is the specific internal energy, and 
q is the energy deposition rate per unit mass. These equations 
assume that the mass of the wind is negligible. Although not 
readily apparent in the form above, eqs. (l)-(3) exhibit a crit- 
ical point when v equals the local speed of sound. In order to 
make the solution to this system tractable and the critical point 
manifest we recast the equations as 
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In the above expressions, C v is the specific heat at constant vol- 
ume, c s is the local adiabatic sound speed, cj is the isothermal 
sound speed, v e = (2GM/V) 1 / 2 , M is the protoneutron star grav- 
itational mass, and 



D = & 



T dP 



e + P dT 



(8) 



Taking the limits v/c <C 1 and c s /c -C 1 we recover the Newto- 
nian wind equations in critical form; 

dv v ( v 2 -4c 2 \ ,. ., 



D q 

dr~ 2r\ c 2 -v 2 J ' Cy~T c 2 -v 2 ' 

dp _ 2p fv 2 -v 2 e /4\ p D q 
dr i 



vC v T c 2 -v 2 



(10) 



and 



dr r C v \ cj-v z J Cyv \ cj-v z 



In this limit, D becomes (T/p) dP/dT\ p . The differential wind 
equations, both Newtonian and general-relativistic, are solved 
in precisely the form above, as a two-point boundary value 
problem, using a relaxation algorithm described in §3. Note 
that direct integration of the continuity equation (eq. 1) yields 
the eigenvalue of the steady-state wind problem, the mass out- 
flow rate M = 4irr 2 pvy. In addition, for q = 0, eqs. (l)-(3) admit 
a second integral of the flow, the Bernoulli integral. We impose 
neither as a mathematical constraint in solving the system. In- 
stead, we use the degree to which each is conserved to gauge 
the accuracy of our converged models. In the Newtonian case, 
the Bernoulli integral can be expressed by 



MA e + 



1 2 P 

-v + — ■ 



GM 



d 3 r'pq = Q(r), 



(12) 



where R v is the coordinate radius of the protoneutron star sur- 
face. In general relativity, with q = 0, 7/1^/— goo is a constant. 
Here, 7 is the Lorenz factor and h is the specific enthalpy. With 
a source term, the differential change in neutrino luminosity is 
given by 



dp 



(13) 



where dp/dr = 4irr 2 pe A and ds 2 = -e 24 'dt 2 + e 2A dr 2 + r 2 dQ, de- 
fines the metric. The total energy deposition rate is then, 

/■OO 

Q = 4tt drr 2 pqe K e 24 >. (14) 

Using these prescriptions for the Bernoulli integral and the 
equation for M and employing modest radial zoning (2000 
points), we typically conserve both to better than 0.1%. 

In our solution to the wind problem, we couple the three wind 
equations in critical form to the differential equation describing 
the evolution of the electron fraction, Y e , due to the charged- 
current electron-type neutrino interactions with free nucleons: 
v e n <-> e~p and D e p <-> e + n. This differential equation does not 
contain a critical point and can be written as 

(vy)-^=X n [T^ n + T e+ „]-X p [T PeP + T e - p l (15) 

where X n and X p are the neutron and proton fraction, respec- 
tively. The T's are the number rates for emission and ab- 
sorption, taken from the approximations of Qian & Woosley 
(1996). The number rate subscripts denote initial-state par- 
ticles. The asymptotic value of the electron fraction (Y") is 
generally determined within ~10 km of the protoneutron star 
surface. Ignoring the details of transport and neutrino decou- 
pling near the neutrinospheres, Y" is determined by both the lu- 
minosity ratio LvjL Ve and the energy ratio (e^) / (e^), where 
(e„) = (El)/(E U ), and E v is the neutrino energy. To rough 
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(but useful) approximation (Qian et al. 1993; Qian and Woosley 
1996), 



Y a ~ 



~ 1 + 



L De ( £Pe )-2A+1.2A 2 /(e Pf ) 



T v , n + T B , P \' ■ L Ve (e,,.)+2A + 1.2A 2 /<e„,) 

(16) 

The energy threshold (A = m n - m p ~ 1 .293 MeV) for the f e 
neutrino absorption process, P e p — » ne + , is manifest in eq. (16). 
Note the difference in sign on the 2A term between numerator 
and denominator. This implies that simply having Lp e jL Vc > 1 
and (sD e )/(Sv e ) > 1 is not sufficient to guarantee Y" < 0.5; the 
magnitudes of (e„ e ) and (e Ve ) are also important. For exam- 
ple, taking UjL Vt = 1.1 and (e Pe }/(e Ue ) = 13MeV/10MeV 
gives Y° ~ 0.52 Hence, if L„ e and L Ve are correlated with 
(sv e ) and (£„,), respectively, then even for constant L^ e jL Ve 
and (ep e ) /{s Ve ), Y" must increase as the total neutrino luminos- 
ity of the protoneutron star decays in time. This phenomenon, 
which we refer to as the threshold effect, is important for the vi- 
ability of an r-process epoch in protoneutron star winds at late 
times and low (e Pf ) and (s Ve ). Other possible effects that might 
bear materially on Y° include the formation of alpha particles 
from free nucleons (the alpha effect) (Fuller and Meyer 1995; 
McLaughlin, Fuller, and Wilson 1996), the differential redshift 
of v e neutrinos versus v e neutrinos, due to the physical separa- 
tion of their respective neutrinospheres (Fuller and Qian 1996), 
and charge conjugation violation (Horowitz and Li 2000). 

Coupled to the wind equations and the equation for Y e evo- 
lution is a simple equation of state (EOS) well-suited to the 
conditions in the neutrino-driven wind (T < 5 MeV, p < 10 13 
g cm" 3 , and 0.0 < Y e < 0.5). Under these conditions, to good 
approximation, free neutrons, protons, and alpha particles may 
be treated as non-relativistic ideal gases. A fully general elec- 
tron/positron equation of state is employed. Photons are also 
included. Past wind studies have approximated electrons and 
positrons as non-degenerate and relativistic. Although r] e (= 
He/T) divided by n may approach ~10 at the protoneutron sur- 
face, it drops steeply with the density so that in the main heating 
region the non-degenerate assumption is justified. In contrast, 
however, it is important to include the non-relativistic character 
of the electrons and positrons. For a broad range of protoneu- 
tron star characteristics, the temperature drops to <~ 0.5 MeV 
within <~ 50- 100 km of the neutron star surface. The dynam- 
ics and asymptotic character of the wind, including mass out- 
flow rate, asymptotic velocity, and composition can be signifi- 
cantly affected by assuming relativistic electrons and positrons 
throughout the wind profile. In addition, the important range 
of matter temperatures for r-process nucleosynthesis occurs for 
T(f) < 0.5 MeV. Sumiyoshi et al. (2000) have found that us- 
ing a general electron/positron EOS can decrease the dynamical 
timescale in the nucleosynthetic region of the wind by as much 
as a factor of two, a potentially important modification when 
considering the viability of the neutrino-driven wind as a can- 
didate site for the r-process. 

3. NUMERICAL TECHNIQUE 

Past numerical studies of steady-state protoneutron star 
winds have at times been hampered by unphysical boundary 
conditions and ill-defined numerics not generally suited to the 
solution of the protoneutron star wind problem. 

The solution to the wind equations constitutes an eigenvalue 
problem. For a given set of protoneutron star characteristics 



and boundary conditions, there exists a unique mass outflow 
rate (M = 4irr 2 pvy) and critical radius (R c ) where the matter ve- 
locity is equal to the local speed of sound (v(R c ) = c s ). Although 
eqs. (5)-(7) are ordinary differential equations, one cannot treat 
the wind as an initial value problem; the critical point necessi- 
tates a two-point boundary value prescription. 

Even though shooting methods determine M precisely, R c 
can remain uncertain (London and Flannery 1982). In addition, 
although Duncan, Shapiro, & Wasserman (1986) made effec- 
tive use of this method, they were forced to employ two types of 
shooting: one for high M solutions and another for low M solu- 
tions. In an effort to circumvent these problems and to construct 
a robust algorithm with flexible boundary conditions, we solve 
both the Newtonian and general-relativistic wind equations as a 
two-point boundary value problem using a relaxation algorithm 
on an adaptive radial mesh, as described in London & Flannery 
(1982) (see also Kippenhahn, Weigert, and Hoffmeister 1968; 
Eggleton 1971; Press et al. 1992). 

The relaxation algorithm involves replacing the differential 
equations with algebraic difference equations at each point on 
the radial grid. Then, given an initial guess for each variable, 
at each point, the solution is obtained by iteration. However, 
because the condition v(R c ) = c s defines the outer boundary and 
this point is not known a priori - such knowledge would effec- 
tively constitute the solution - we follow the procedure of Lon- 
don & Flannery (1982) and introduce a new independent vari- 
able q that labels radial mesh points by integers (q\<q< qu). 
The price paid is three more differential equations: 

dr ip 
dq 4>{r) 1 



dq 



and 



dip 
dq 



= 0. 



(17) 



In this scheme, r becomes a dependent variable, Q(r) is a mesh 
spacing function (e.g., Q(r) = logr), ip is an intermediate vari- 
able, and <fr(r) is proportional to the density of mesh points. 
The system is solved on the mesh of q values. Hence, the outer 
radial coordinate, at which the v(R c ) = c s boundary condition is 
to obtain, adjusts in a Newton-Raphson sense to simultaneously 
satisfy all boundary conditions. Typically, when we begin a cal- 
culation, we start with an initial guess that extends to a radius 
with Mach number of ~ 0.9. With each iteration, R c is adjusted 
so that the Mach number goes to 1 . 

We now have three wind equations for p(q), T(q), and v(q), 
three extra differential equations for the mesh algorithm, and 
the equation for Y e (q). Seven boundary conditions must then be 
imposed to close the system. The boundary conditions on the 
first two of eqs. (17) are simply r{q{) = R v and Q(q\) = \og(R u ) 
(for log spacing), where R v is the protoneutron star radius. Two 
more boundary conditions obtain at the critical point. Taking 
qti as the outer mesh point, we have that v(qrf) = c s and we have 
that the numerator of any of the differential equations for p(r), 
T(r), or v(r) must simultaneously be zero to ensure continuity 
of the solution through the sonic point. Three additional bound- 
ary conditions are required to specify the problem completely. 
Although simply setting Y e (q\), p(q\), and T(q\) at R„ is suf- 
ficient, we do not use this prescription. Instead, we assume 
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that the radius of neutrino decoupling coincides with the coor- 
dinate radius of the protoneutron star surface. In addition, we 
assume that this neutrinosphere (/?„) is the same for all neutrino 
species: electron (v e ), anti-electron (D e ), and mu and tau neutri- 
nos (collectively, i^'s). Indeed, as the protoneutron star cools 
we expect these neutrinospheres to be separated by just tenths 
of kilometers. As the v e neutrinos have the largest net opacity 
(n Ve ) of any species, we set an integral boundary condition on 
their optical depth (r Ve ): 



T„ e (R v ) 



n Ve pdr = -. 



(18) 



Included in the v e opacity are contributions from scattering with 
free nucleons, scattering on electron/positron pairs, v e n — > pe~, 
and alpha scattering. As a second boundary condition, we as- 
sume that the net neutrino heating balances the net cooling at 
R v . That is, q{R v ) = 0. Finally, our third boundary condition 
assumes that the charged-current processes are in equilibrium 
at the protoneutron star surface. Explicitly, 

^ =0. (19) 
dr ** 

For a given protoneutron star mass and radius and a set of av- 
erage neutrino energies and luminosities, we start with a guess 
for the mass density at the surface of the protoneutron star. We 
use a two-dimensional Newton-Raphson algorithm to simulta- 
neously satisfy the conditions on q and dY e /dr at R v . This de- 
termines T(R V ) and Y e (R v ) for the first step. At each subsequent 
iterative step, the relaxation algorithm attempts to satisfy the 
integral boundary condition on r„ c . Effectively, this procedure 
results in a new T(R U ) and Y e (R u ). In this way, we satisfy all 
boundary conditions simultaneously. Given a good initial guess 
for the solution (i.e., maximum deviations from convergence in 
any variable of < 20%) we obtain a solution in just 5-10 iter- 
ations. Once the profile for /?„ < r < R c is obtained, we use 
l'Hospital's rule to bridge R c and then integrate to larger radii 
as an initial value problem using a fourth-order Runga-Kutta 
scheme. Successfully converged models then serve as an ini- 
tial guess for the next protoneutron star model with adjacent 
characteristics (i.e., in mass, radius, or neutrino spectral char- 
acteristics). 

3.1. Tests of the Code 

We do not impose dM/dr = or the Bernoulli integral as 
mathematical constraints on the system. Instead, we use the de- 
gree to which these conditions obtain to gauge the precision of 
our method. Typically, both are conserved to better than 0.1% 
in both general-relativistic and Newtonian calculations. In ad- 
dition, as we increase the number of mesh points, the error in 
both of these quantities decreases significantly. 

One may argue that eq. (19) need not hold generally. In fact, 
we adjusted the code to accept a fixed Y e (R v ) boundary condi- 
tion to make sure this has no effect on the asymptotic character 
of the wind. We find that the number rates are large enough 
at the surface that Y e is forced from Y e (R v ) to a value such that 
dY e /dr <~ in the first radial zone with no appreciable effect on 
any aspect of the wind. For simplicity, then, we have chosen to 
enforce eq. (19). 

Similarly, in a fully dynamical calculation, one would not 
expect q{R v ) = 0, generally. While this may certainly be true, 



in our solution to the steady-state wind we encounter numeri- 
cal instabilities that preclude solution of the equations with fi- 
nite q{R v ). It is difficult to estimate the importance of such an 
assumption without employing the full machinery of radiation 
hydrodynamics. 



4. THE NEUTRINO HEATING FUNCTION 

The neutrino energy deposition rate (q) is a sum of contri- 
butions from the charged-current v e and v e neutrino absorp- 
tion processes, neutrino-electron/positron scattering, neutrino- 
nucleon scattering, and the process vv <-> e + e~. We describe 
each in turn. 



4.1. The Charmed-Current Processes 

At the entropies encountered in supernovae (< 40), the 
charged-current or beta processes (y e n <-» e~p and v e p <-> e + n) 
dominate the opacity and energy exchange for the electron and 
anti-electron neutrinos. In the protoneutron star wind context, 
at much higher entropies, we expect these processes to compete 
with neutrino-electron/positron scattering in the net energy de- 
position. Ignoring final-state blocking and assuming relativistic 
electrons and positrons, the charged-current cooling function 
can be written as 



C cc ~2.0x 10 18 T 6 



Xp F 5 (0) +Xn F 5 (0) 



erg g 1 s 1 , 



(20) 



where 



F„(y) = 



o e*-y+i 



-dx 



are the Fermi integrals, X„ and X p are the neutron and proton 
fractions, respectively, T is in MeV, and t] e = /i e /T. The heat- 
ing rate due to neutrino captures on free nucleons can be written 

as 



s" 1 



ff cc ~ 9.3xl0 18 ^ 6 2 ergg 

x [X n Ll](el e )+X p Ll](el e )] $ 6 H(r), 



(21) 



where R^ is the neutrinosphere radius in units of 10 6 cm and 
both L v and (e„) are defined at R v . We separate L v and (e„) 
in this heating rate and those below. Although L v and (e v ) are 
correlated, L v and the tail of the neutrino energy distribution 
are generally not. We retain L v and (e„) separately so we have 
the freedom to change them independently. With eqs. (20) and 
(21), the net energy deposition due to the charged-current pro- 
cesses is then q cc = H cc -C cc . In eq. (21), 



(el) ^ 



J d[i J de u elf„ ■ J d\i J 



(22) 



where /i(= cos 9) is the cosine of the zenith angle and /„ is the 
neutrino distribution function. $ accounts for the gravitational 
redshift of neutrinos from the protoneutron star surface and is 
given by 

\ l-lGM/rc* ) ■ ( } 

Also present in eq. (21) is the spherical dilution function, S(r), 
which describes the radial dependence of the neutrino energy 
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and number densities. In the vacuum approximation, assuming 
a sharp neutrinosphere, 



E(r)=l-Jl-(R,/r) 2 /^ 



(24) 



and is related to the flux factor (/i)(= F v /J v ) by (fi) = 
/? 2 /2$ 2 S(r)r 2 , where F v and J v are the neutrino flux and en- 
ergy density, respectively. In eq. (24), the factor $ 2 accounts 
for the bending of null geodesies in a curved spacetime (Cardall 
and Fuller 1997; Salmonson and Wilson 1999). This effectively 
increases the neutrino number density a given mass element 
sees at any radius, thus augmenting the energy deposition. In 
contrast, since the heating rate for any neutrino interaction is 
proportional to positive powers of L v and {ej), the gravitational 
redshift terms which modify these quantities can only decrease 
the energy deposition rate at a given radius. As pointed out by 
Cardall & Fuller (1997), the augmentation of q by the bending 
of neutrino trajectories and the degradation of q by the gravita- 
tional redshift compete as M/R u increases, but the latter domi- 
nates. 

Qian & Woosley (1996), Otsuki et al. (2000), and Wanajo 
et al. (2000) have employed eq. (24) or its Newtonian analog 
in their studies of neutrino-driven winds. In an effort to ad- 
dress neutrino decoupling more fully, we compared our wind 
solutions using this spherical dilution factor to those obtained 
with an effective S(r) derived from the Monte-Carlo transport 
results of Janka (1991) who connected (/i) with the density 
gradient at 10 10 g cm" 3 and the curvature of the opacity pro- 
file at the radius of decoupling. The difference between the 
radial dependence of (/i) using this approach and the effective 
flux factor obtained in a vacuum approximation is significant, 
as would be expected, only when the atmosphere is sufficiently 
extended. That is, as the density gradient just exterior to the 
neutrinosphere goes to infinity, so too should the Monte-Carlo 
results of Janka (1991) approach the vacuum approximation. 
Hence, significant deviations from the vacuum approximation 
only present themselves when the radius of the neutron star is 
large (~ 20-40 km) and/or the total neutrino luminosity is very 
high (i.e., the temperature at R v is large). As the protoneutron 
star cools and the luminosity decreases for a given R v , the den- 
sity gradient becomes steeper (see Fig. 2 below), thus making 
the vacuum approximation more appropriate. Over the range of 
models presented here, we have examined the effects of using 
the Janka (1991) (/z) in order to characterize the spherical di- 
lution of neutrinos for all of the relevant heating processes and 
find them negligible, particularly for the compact protoneutron 
star wind models most likely to be important for r-process nu- 
cleosynthesis. 

We emphasize that this simple comparison is not a complete 
analysis of the transport effects that might be important in de- 
termining q(r). While this result suggests the vacuum approx- 
imation is appropriate in our scheme for handling the neutri- 
nospheres and boundary conditions, it says nothing about actual 
transport effects in the decoupling region. 



4.2. Neutrino-Electron/Positron Scattering 

At high entropies, electron-positron pairs are produced in 
abundance. Therefore, neutrino-pair scattering is expected to 
contribute substantially to the total energy deposition in the 
protoneutron star wind. The energy transfer associated with a 



single neutrino-electron or positron scattering event is well ap- 
proximated by Ui ~ (e Vi —AT)/2, where iv is the energy trans- 
fer, T is the matter temperature, and e Vt is the neutrino energy 
of species ;' (Bahcall 1964). We have confirmed this approxi- 
mation using the thermalization code described in Thompson, 
Burrows, & Horvath (2000), which employs the fully relativis- 
tic structure function formalism of Reddy et al. (1998). The net 
heating rate due to the interaction of the neutrino fluid with the 
pair plasma can be approximated by q ~ cn e n Vi (a Vte oj), where 
n e and n Vi are the number density of electrons and neutrinos, 
respectively, and a v<e is the cross section for scattering. We 
obtain 

/dn e dn v . 
w c o Vie ——de e - — -de v , , (25) 
de e de Vi 

where a Vie (= T e Vi ) is the cross-section for neutrino scattering 
on relativistic, non-degenerate electrons (Tubbs and Schramm 
1975). Ki = cr A,'/2m 2 is a neutrino species dependent constant, 
where m e is the mass of the electron in MeV, a Q ~ 1 .71 x 10" 44 
cm 2 , and 

1 



Ai = (c v +c A y 



+ ^(c v -ca) 



(26) 



cy and c A are the vector and axial-vector coupling constants for 
a given neutrino species. We find that the energy deposition rate 
can be expressed as 

u 

erg g 



c ( T F 2 (r le ) 



K 

2 {£ 



(He) 3 TT 2 

'F 3 (r)„) J 



47rr 2 c(£„)(/i) 
F 2 (Vv) 



$ 2 



F}{T}v} 



AT 



FliVe) 



■(27) 



We drop the subscript i here for brevity. Note that (ii)(= 
R 2 /2$ 2 S(r)r 2 ) contains two powers of the redshift. r\ v is an ef- 
fective degeneracy parameter obtained from fitting the neutrino 
distribution function with a Fermi-Dirac distribution with an 
appropriate T v and r\ v (Janka and Hillebrandt 1989). Although 
we retain the general form here, in the wind calculations pre- 
sented below we take r\ v = for all neutrino species. In eq. (27), 
the first term in parentheses is the number density of electrons, 
n e . Alternatively, with the replacement j] e — ► — r\ e , it is the num- 
ber density of positrons. The next term is the number density 
of neutrinos, a function of radius, which depends also on the 
flux factor, (fi). The term in square brackets is the appropri- 
ately averaged product of the cross section for scattering and 
the energy transfer per scattering, (a v , e u>) . We retain the form 
q = cn e n Vj (a Uie u>) here for clarity. Note that, for i] u = r\ e = 0, 
(e„) equals 3.1571, and the second quantity in parentheses in 
eq. (27) is simply (T Vl - T); net heating occurs if T v > T(r) at 
any r. In order to obtain the contribution to the net heating from 
neutrino-positron scattering, in addition to the change in r] e , one 
must also make appropriate changes to A,. 

4.3. Neutrino-Nucleon Scattering 

The energy transfer associated with a single neutrino- 
nucleon scattering event is much smaller than that for neutrino- 
electron scattering. The cross section, however, is much larger. 
Using our neutrino thermalization code, which solves the Boltz- 
mann equation in an isotropic, homogeneous thermal bath of 
scatterers, including the full collision term, with Pauli blocking 
and explicit coupling between all energy bins (Thompson, Bur- 
rows, and Horvath 2000), we confirm that the average energy 
transfer for neutrino-nucleon scattering, in non-degenerate nu- 
clear matter, for neutrino energies below ^40 MeV is well ap- 
proximated by w ~ e u (e u -6T)/m N , where m N is the nucleon 
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mass in MeV (Tubbs 1979). We derive the heating rate for 
neutrino -neutron scattering as 



p 4TTr 2 c(e v )(p) 

{£„) F 6 (r)„) 

m„ F 3 (r] v ) 



4> 4 



erg g 



1 s" 1 



F 2 (Vv) 6T F5(v e ) 



FliVv) " F 6 (jie) 

(28) 

k = a /(l6m^)(l + 3g\) for neutron scattering and n = 
<7 /(4m 2 )[4 sin 4 8 w -2 sin 2 9 w + {\ + 3g%) /4] for neutrino-proton 



scattering, where sin 2 8 W ■ 
vector coupling constant. 



: 0.231 andg A (— -1-26) is the axial- 



We were motivated to consider nucleon-nucleon 
bremsstrahlung, plasma, and photo-neutrino processes by the 
sensitivity of the dynamical timescale and asymptotic entropy 
to changes in the energy deposition profile (see §8.3 and Qian 
and Woosley 1996). We found that none of these processes 
contributed significantly. Qian and Fuller (1995a,b) have ad- 
dressed the possibility that neutrino oscillations may effectively 
provide an extra energy source, at larger radii, beyond the point 
where the mass outflow rate is determined. In §8.3, we include 
an ad hoc energy source to test the sensitivity of our results 
to changes in q(r), but do not address the issue of neutrino 
oscillations directly. 



4.4. ViVi <-> e + e 

Cooling due to e + e~ annihilation, assuming relativistic elec- 
trons and positrons, and ignoring Pauli blocking in the final 
state, can be written as 



~ 1.4 x 10 17 — f(T) e ) ergs g 1 s 



where 



f(Ve) = 



F 4 (T] e )F3(-T] e ) + F4(-T]e)F3(7] e ) 
2F 4 (0)F 3 (0) 



(29) 



(30) 



P8 is the mass density in units of 10 8 g cm" 3 and T is in MeV. 
We employed eq. (29) with and without the r\ e dependence and 
compared it to the results using the fit given in Itoh et al. (1996). 
Such a modification amounted to no more than a 1-2% change 
in the entropy, dynamical timescale, total energy deposition, or 
mass outflow rate. 

Much more uncertain than eq. (29) is the heating due to 
vv — > e + e~. The spherical dilution of this process from the neu- 
trinospheres is complicated. For a given flux, the local energy 
density depends sensitively on (/i), which can only be prop- 
erly treated by solving the full transport problem. We save such 
an investigation for a later work. Instead, we compare two ap- 
proaches. The first is based on the vacuum spherical dilution 
approximation, written simply as (Qian and Woosley 1996) 

H~1.6xl0 19 ^S- ^ergg" 1 s" 1 



L^Ll%e, e ) + {e Vc )) + -{Liy{e^) 



(31) 



where = (1 -x)\x 2 +4x+5), x = (1 -{R v /r? /$> 2 ) l l 2 , R v6 
is the neutrinosphere radius in units of 10 6 cm, and is the 
neutrino luminosity in units of 10 51 erg s" 1 . The redshift term, 
$, appearing in x, accounts for the amplification of this process 
due to the bending of null geodesies in general relativity as in 
eq. (24) (Salmonson and Wilson 1999). We compared this ap- 
proximation to the heating rate obtained by Janka (1991). Be- 
cause of the extreme nature of the density gradient just exterior 
to the protoneutron star, we find that the vacuum approximation 
adequately characterizes the energy deposition. In addition, the 
fact that there is no obvious way to include the general relativis- 
tic effects in the parameterization of (/i) by Janka (1991) led 
us to employ eq. (31) in the wind models we present here. 

4.5. Other Possible Neutrino Processes 



5. RESULTS: FIDUCIAL MODELS 

With a robust and efficient means by which to solve the 
wind problem, coupled with physical boundary conditions and 
a well-motivated neutrino heating algorithm we can now survey 
the protoneutron star wind parameter space. 

Supernova and protoneutron star calculations, coupled with 
observations of neutron star binary systems and our knowledge 
of the high density nuclear equation of state, place useful limits 
on the parameter space that protoneutron stars actually inhabit. 
Particularly, we are interested in a range of protoneutron star 
masses from 1.2 to 2.0 M©, total neutrino luminosities from 
4 x 10 52 to 1 x 10 51 erg s" 1 , average neutrino energies as high 
as 35 MeV for neutrinos and 15-20 MeV for v e and v e neu- 
trinos with (e Ve ) < (sp e ) < {e„ ), and a range of neutron star 
radii from 20 km to perhaps 9 km. 

Our main goal in this section is to map the possible pro- 
toneutron star wind parameter space, taking as input the phys- 
ical ranges specified above. In what follows, we present the 
basic structure of the wind and identify some of the systemat- 
ics we can use to assess this site as a candidate for r-process 
nucleosynthesis. Particularly, we include a discussion of the 
dynamical timescale, the electron fraction, and the asymptotic 
entropy. In identifying some of the approximate power law re- 
lations which characterize the wind, we refer to low and high 
total neutrino luminosities, with low denoting < 10 52 erg s -1 
and high meaning > 10 52 erg s" 1 . We write these power laws 
in terms of L v and (e v ). These stand for a representative lumi- 
nosity and a representative average energy, respectively. That 
is, we take L D jL Ve , L P JL^, (s De )/{e Ve ), and (£*„)/(£„„) as 
constant. 

We employ the notation L 5 J to refer to luminosities in units of 
10 51 erg s" 1 . All neutrino luminosities quoted throughout this 
paper are local quantities at the neutron star surface. The lumi- 
nosity at infinity can be obtained from the luminosities quoted 
here by multiplying the local luminosity by two powers of the 
gravitational redshift (eq. 23). For example, taking M = 1 .4 M Q 
and R v = 10 km, $ 2 ~ 0.58. In order to keep track of our mod- 
els, we use the v e neutrino luminosity (Lp e ) to label each indi- 
vidual model. The total luminosity can then be obtained from 
the ratios L D jL Ve and L Pe /L^. Finally, in our expressions for 
q we included r\ Vj terms for completeness. In what follows, we 
take r\ Ve = r\ Ve = r\ v =0. For an assessment of this assumption, 
see Janka & Hillebrandt (1989) and Myra & Burrows (1990). 
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5.1. The Structure of Neutrino-Driven Winds 

Figures 1, 2, and 3 show the velocity, mass density, and 
temperature as a function of radius for eight different neu- 
trino luminosities. For these models, R u is 10 km and the 
neutron star gravitational mass is 1.4 M Q . The average neu- 
trino energies were set at (s Ue ) = 11 MeV, (e Pf ) = 14 MeV, and 
(e v ) = 23 MeV for the highest luminosity. For each subse- 
quent luminosity, the average energies were decreased accord- 
ing to oc L l J A . The luminosities were held in the ratios 
LvjL Ue = 1.3 and L„ e jL v ^ = 1 .4. In these figures, critical points 
(v = c s ) are shown as dots. Table 1 lists the global properties of 
each of these models, including the integrated energy deposi- 
tion rate, Q (eq. 14), and the mechanical luminosity or hydro- 
dynamical power, P mec h = Mv 2 a /2, where v a is the asymptotic 
velocity. In addition, we include the asymptotic entropy s a and 
dynamical timescale defined as the e-folding time of the density 
at T = 0.5 MeV: 



Figure 1 shows clearly that even relatively close to the neu- 
tron star (r <~ 400 km) the flow is not homologous (i.e., v is 
not proportional to r). At much smaller radii, however, in the 
heating region (r < 60 km, compare with Fig. 4) the flow is ho- 
mologous with v oc L v r. We have found that simple parameter- 
izations of v(r), particularly for low neutrino luminosities, are 
not straightforward. Although not readily apparent, the asymp- 
totic velocity (v a ) is a power law in L v . For high L}°\ v a oc L? v 3 . 
For low L u , the index is ~ 0.46. The critical radius (R c ) also in- 
creases as a power of luminosity for low L v . For M = 1 .4 M©, 
the index is - 0.95, while for M = 2.0 M R c oc L"° 85 . Note 
that if the protoneutron star cools as L v oc f "° 9 then we should 
expect R c to grow linearly with time, implying that at late times 
the velocity of the critical point away from the neutron star is 
approximately constant. 

The stiff nature of the wind equations is manifest in Fig. 2. 
The inset shows that p drops by more than four orders of mag- 
nitude in just 1-3 km, before the neutrino heating rate reaches 
a maximum. Figure 4 shows the total specific heating rate for 
each of the eight fiducial wind models. Note that the maxi- 
mum in q occurs very close to the protoneutron star surface 
and that the position of this maximum is not a function of L v . 
At r <~ 12 km the density gradient changes dramatically, as the 
slow moving material in this inner atmosphere accelerates to 
infinity. Concomitant with this change in p(r) and the peak in 
q is a rise in T and most of the entropy production. Figure 5 
shows the entropy as a function of radius. Note that the en- 
tropy quickly comes to within a few percent of its asymptotic 
value (s a ). We find that s a is proportional to L~ 24 for all lu- 
minosities, for 1.4 M . For M = 2.0 M Q , s a is proportional to 
L~ 25 . As L v decreases, the gradients in p, T, and the entropy 
near /?„ become larger. Inspection of Fig. 3 shows that the 0.1 
MeV < T < 0.5 MeV region where the r-process (and preced- 
ing a-process) might take place, lies between 50 km and 500 
km for L%\ =8.0 and between 25 km and 150 km for L%] = 1 .0. 
Importantly, Fig. 1 shows that the wind is still accelerating sig- 
nificantly at these radii. This calls into question the assumption 
of a constant velocity outflow used by many r-process modelers 
(e.g., Meyer and Brown 1997). 

Figure 6 shows a breakdown of the total energy deposition 
into its separate components for the L 5 ^ =8.0 fiducial model. 



All of the heating processes conspire to give q(R u ) = 0, in accor- 
dance with our boundary condition. At the surface, the charged- 
current heating rate (q cc ) provides net cooling, which balances 
heating from neutrino-nucleon scattering (q v n)- q C c generally 
dominates the other heating processes in the models we con- 
sider here. Note that heating due to neutrino-electron/positron 
scattering (q ve ) does not contribute as significantly as is indi- 
cated by Qian & Woosley (1996) and Otsuki et al. (2000), 
because of their simplifying assumption that the energy trans- 
fer per scattering is <~ e v /2 instead of <~ (e„-4r)/2 (§4.2). 
Comparing Fig. 6 with Fig. 7, we can see that the drop in 
q cc at r ~ 35 km is due to the formation of a particles. As 
the luminosity decreases this transition region moves in, so 
that for Lp' = 0.60, q cc — ► near r <~ 22 km. As the ratio of 
MjR v increases, the heating rate due to — > e + e~ {q V v) in- 
creases substantially relative to both q ve and q cc . In fact, for 
MjR v = 2.0M Q /10 km, the peak in q vV is actually slightly 
higher than that for q ve . However, q cc still dominates heating at 
the peak in q\ ot by a factor of <~ 1 .5. Note that an increase in L„ M 
or (e„ ) increases the importance of q ve and q V v relative to q cc . 
In particular, this increases the energy deposition at larger radii. 
We explore some of these effects on the asymptotic character of 
the wind in §8.3. 

Qian & Woosley (1996) found that M oc L 2 - 5 AT 2 for Newto- 
nian gravity. We find that the index of this power law in L v is 
slightly decreased in general relativity to 2.4-2.5, but that the 
dependence on M is stronger; M is approximately proportional 
to M~ 3A for high luminosities. The index decreases to -2.7 for 
low luminosities. For the full range of luminosities and masses 
considered here, the mass outflow rate in our general-relativistic 
calculations is approximately a factor of three smaller than in 
our Newtonian wind models. 

Computing the volumetric integral of q as in eq. (14) yields 
the net energy deposition rate, Q (see Table 1). For the heat- 
ing function we employ in this paper, Q is roughly proportional 
to L 2,4 for all protoneutron star masses and luminosities. The 
small variations in this power law index compliment the vari- 
ation in M with L v so that the ratio Q/M varies by less than 
2% over the whole range of L v for each mass. In all cases 
Q/M ~ GM/R V to within 10%. Whereas Qian & Woosley 
(1996) found that Q oc M" 1 , we find a more stiff dependence 
on mass; for high L v , Q is proportional to M~ 19 for 1.4 M Q 
< M < 1.6 M Q . In addition, we find that their analytic expres- 
sion for Q consistently underestimates the net energy deposition 
in our Newtonian models by as much as ~ 50%. 

The wind mechanical power P mec h is proportional to L 3 U 4 for 
lowL„ andM= 1.4M Q . This index decreases to3.2forM = 2.0 
Mq. At high luminosity for all masses, P mec h oc L]; 1 . Although 
the asymptotic velocity is set by the escape velocity from the 
protoneutron star and therefore increases with increasing M, it 
does not increase as a power law in general relativity. In fact, 
it increases more rapidly. Therefore, even though M is approx- 
imately proportional to M" 2 5 at low luminosities, the increase 
in v a as M gets large forces P mec h to increase. 

The efficiency of energy deposition, Q/I}° 1 , ranges from 10~ 3 
to 10" 5 as L v decreases. The efficiency of conversion of neu- 
trino energy to hydrodynamical power, Pmech/i 1 "', ranges from 
only 2 x 10" 5 to less than 6 x 10" 8 for the models considered 
here. The efficiency for the conversion of net energy deposition 
to hydrodynamical power, Pmech/2, decreases with luminosity 
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as I} V A for 1.4 M Q with P mech /Q ~ 2.4 x 1(T 3 for = 0.6. 
This means that almost all of the net energy deposition goes 
into overcoming the gravitational potential. The excess energy, 
manifest at infinity as the mechanical power, is very small in 
comparison with Q. These quantities may be potentially im- 
portant if the wind is to emerge and escape to infinity in the 
expanding supernova envelope. 

5.2. The Effects of General Relativity 

Over the range of masses presented here, we find significant 
enhancements in the entropy per baryon using the full general- 
relativistic framework. Over a broad range of luminosities for 
the 1 .4 M Q protoneutron star we find that s a is 25-30 units less 
in our Newtonian calculations than in our analogous general- 
relativistic calculations. Typical reductions in M and Q are 
of order a factor of three and two, respectively. These differ- 
ences were anticipated by Qian & Woosley (1996) and Cardall 
& Fuller (1997) and more recently realized in the wind cal- 
culations of Otsuki et al. (2000). As the latter showed, the 
general-relativistic effects on s a are much more the result of us- 
ing the general-relativistic hydrodynamic formulation than of 
incorporating the general-relativistic corrections to q expressed 
in eqs. (23) and (24). The inclusion of general relativity in 
the hydrodynamics makes the structure of the protoneutron star 
more compact than in a Newtonian description. This makes 
the temperature and density gradients steeper just exterior to 
R u and particularly in the heating region. Although dT jdr 
decreases rapidly, dp jdr drops much faster. This effectively 
increases the specific energy deposition per unit mass and the 
entropy is enhanced significantly. 

For comparison, we calculated several wind models with 
eqs. (5)-(7) that included the enhancement to q due to the bend- 
ing of null geodesies, but did not include any redshift factors on 
L„ or (e v ). For the 1.4 M model with L 5 ^ = 8.0 (see Table 1) 
this change nearly doubled both Q and M, while decreasing r p 
from 3.68 to 1 .95 ms. Although more heating occurred at larger 
radii, the peak in q also increased so as to offset any potential 
gain in s a . The net enhancement was just 6 units in entropy. 
Conversely, keeping all the redshift terms and eliminating <f> 
from E(r) in eq. (24) increased t p by more than 30% and left s a 
virtually unchanged. 

Because the bending of null geodesies increases the net en- 
ergy deposition close to the neutron star and the redshift terms 
act to decrease energy deposition over the whole profile, it is 
clear from this comparison that only an increase in q(r) that 
does not significantly increase Q can have large effects on s a 
(see §8.3 and Qian and Woosley 1996). 

5.3. The Electron Fraction 

Figure 7 shows the evolution of the neutron (X n ), proton (X p ), 
alpha (X a ), and electron fraction (Y e ) as a function of radius for 
the Ll] = 8.0 fiducial model (Table 1 and Figs. 1-6). The elec- 
tron fraction profile is computed using eq. (15), solved simul- 
taneously with the wind equations. In computing the wind so- 
lutions we assume nuclear statistical equilibrium between free 
nucleons and alpha particles. Y e comes to within a few percent 
of its asymptotic value (Y") in just the first five kilometers. This 
quick evolution is due primarily to the low matter velocities in 
the inner region. For L^ e /L Ve =8/6.15 and (s Pe )/{e Ue ) = 14/11, 



eq. (16) predicts that Y" ~ 0.478. Solving the differential equa- 
tion we find remarkable agreement: Y e ~ 0.477 at r = 20 km. 
Just beyond this, for r > 30 km, free nucleons form a particles 
and Y e rises slightly as X a increases until Y e = 0.485 at r ~ 150 
km. This is the a effect, whose import in this context was first 
noted by Fuller & Meyer (1995) and McLaughlin, Fuller, & 
Wilson (1996). The magnitude here is only of order 1%. For 
Y" < 0.5, we find generally that the magnitude of the a effect 
increases as the luminosity decreases for a given R v . In addi- 
tion, for models in which Ljj e jL Ve and (sD e )/(s Ve ) are larger, 
and, hence, Y" is naively lower (a la eq. 16), the magnitude of 
the a effect is also enhanced. However, for a reasonable range 
of L„jL Ve and (s^) j '{e Ve ) as well as total neutrino luminosi- 
ties, the a effect never increases Y" by more than ^10%. That 
is, if Y e = 0.40 before a particle formation, we find that the a 
effect increases Y e to no more than approximately 0.44. 

Since in these fiducial models we decrease the average 
neutrino energies with luminosity, the threshold effect in the 
charged-current reactions, manifest in eq. (16) by the neu- 
tron/proton mass difference (A), becomes important. Despite 
the fact that LyjL Ve and (sa e )/(s Ve ) are maintained as above, 
Y" eventually becomes greater than 0.5. As Y" becomes much 
larger (<~ 0.52), it experiences what might also be termed an a 
effect: because an a particle has Y e = Z/A = 0.50, the onset of 
X a formation decreases Y". 

5.4. The Dynamical Timescale 

In discussing our results, in order to make an apposite com- 
parison with previous studies, we quote t p (eq. 32) at T = 0.5 
MeV. However, using such a scale to characterize the nature of 
the resulting nucleosynthesis is suspect. Figure 8 shows r p as 
a function of radius for the wind models presented in Figs. 1- 
5. The dots on each line of constant neutrino luminosity mark 
the range 0.5 MeV > T{r) > 0.1 MeV, the temperature range 
relevant for neutron-capture nucleosynthesis. Although the r- 
process may continue at temperatures well below 0.1 MeV, we 
include these dots to guide the eye. The dynamical timescale 
of the wind, or the expansion timescale, has been defined in 
several different ways by many researchers. No clear consen- 
sus exists. Cardall & Fuller (1997) defined their dynamical 
timescale as the e-folding time of the temperature at T = 0.5 
MeV. Hoffman, Woosley, & Qian (1997) defined their dynam- 
ical timescale in the same way, but then used an expansion 
timescale (1.28 times the dynamical time) to discuss their re- 
sults. Qian & Woosley (1996) and Freiburghaus et al. (1999) 
used the ratio r/v at T = 0.5 MeV to characterize the expansion. 
Finally, Meyer & Brown (1997) connect the e-folding time of 
the density with their expansion timescale, r/2v, by using the 
equation for M and dropping the acceleration term. With these 
assumptions, they obtain 

t„ mb 00 = to(i + ^), (33) 

where / is the time on a Lagrangean mass element in the flow. 
Figure 8 shows clearly that any simple parameterization of the 
dynamical time, using any definition, is an oversimplification. 
For Lp' = 8, r p increases by a factor of three over this range of 
temperatures. Equation (33) captures the increase of r p with ra- 
dius for high luminosities, but overestimates the slope by about 
a factor of two. At low luminosities, of course, eq. (33) does 
not capture r p (f) at all. At these luminosities, the dynamical 
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timescale actually decreases over this range of temperatures. 
This arises because a region of positive curvature in the v(r) 
profile of Fig. 1 develops between 40 km and 200 km at low 
luminosities. This observation simply underscores the danger 
in considering a single dynamical time that is meant to charac- 
terize an actual wind profile. 

Figure 9 shows tracks of constant mass in the plane of s a ver- 
sus t p for luminosities from L 5 ^ = 8.0 to L 5 ^ = 0.70, for R v = 10 
km. Although the indices vary slightly for different masses, 
s a oc L~ 25 and r p oc L~ 1A , so that s a is approximately propor- 
tional to Tp 2 . One can imagine that these curves represent evo- 
lutionary cooling tracks in time in the space of s a and t p for 
constant R v and M. Simple extrapolation of this power law al- 
lows one to estimate at what t p a given entropy might obtain, 
for a given mass and neutron star radius. For example, the 1.4 
M trajectory in Fig. 9 will not reach s a ~ 200 until r p ~ 0.4 
seconds. The corresponding neutrino luminosity at this point 
is L 5 J e ~ 0.22. Knowing that M oc L 2 ; 4 allows one to estimate 
the mass outflow rate as M ~ 1 .8 x 10" 8 M Q s" 1 . These simple 
power laws and a knowledge of how L v might behave in time 
allow us to put powerful constraints on the likely wind epoch 
of r-process nucleosynthesis, as we demonstrate in §6. 

Wanajo et al. (2000) found that their dynamical timescale 
saturated at high neutrino luminosities near 2 - 3 x 10 52 erg s" 1 . 
This conclusion is an artifact of their definition for the dynam- 
ical timescale and their numerical approach to the wind prob- 
lem. Our solution shows that for constant average neutrino en- 
ergies, even up to L l ° l ~ 8 x 10 52 erg s" 1 , t p continues to de- 
crease, roughly as t p oc L~ 0J5 . The entropy also decreases as 
the luminosity increases; s a oc L~ 0A5 . Since these power laws 
are for constant average neutrino energy, we deduce from the 
fiducial models that t p oc L? 15 (e^)' 2 - 6 and s a oc L? AS {e v )~° A . 
These are to be compared with the analytic results of Qian & 
Woosley (1996) who did not include general-relativistic ef- 
fects: Td yn oc L~ l (e u )~ 2 and s a oc L„ 1/ ' 6 (£ ly )~ 1 < /3 . In addition, al- 
though we find that r p decreases as R v decreases, for constant 
protoneutron star mass, neutrino luminosity, and average neu- 
trino energy, using either the Newtonian or general-relativistic 
wind equations, we find that r p is generally 10%— 15% shorter 
in the Newtonian case. This result is owed in part to the fact 
that the increase in q due to the bending of neutrino trajectories 
is insufficient to counter the decrease in q caused by gravita- 
tional redshift of the neutrino luminosity and energy (Cardall 
& Fuller 1997). (See Table 1 for a comparison between the 
Newtonian and general-relativistic calculations for the 1 .4 M 
and L|J =8.0 fiducial model.) 

5.5. The Steady-State Approximation 

To conclude this section, we include a few words about the 
degree to which the protoneutron star wind can be considered 
quasi-stationary. There are several timescales of importance. 
The first is n, the timescale for decay of the neutrino lumi- 
nosity, set by the power-law index S in the relation L v oc t~ s : 
i~d = t /S. The second is the time t„, for matter to move from 
R v to the critical point R c , where it loses sonic contact with 
the rest of the flow. The third relevant timescale is the sound 
crossing time, t s , between R u and R c . t, varies from just <~ 10 
ms to more than 250 ms over the range of L„s presented here. 
For low luminosities, t s oc L~ 1A . r m is proportional to L~ 2 8 for 
all luminosities and varies from 4 seconds to more than 5000 



seconds for the same range of L„s. However, these numbers 
for r m are quite deceiving. In our models, due to the exponen- 
tial density gradient just exterior to R u , the matter is effectively 
trapped for r < 12 km. In fact, for the lowest luminosity cases 
presented here, the matter velocity for r < 10.5 km can be of 
order 10 cm s -1 . In effect, then, the region shown in the inset in 
Fig. 2 is an atmosphere in hydrostatic equilibrium from which 
the wind emerges. From this region the matter would escape on 
timescales much longer than the total protoneutron star cool- 
ing time. If, instead, we redefine r m as the time necessary for 
a Lagrangean mass element to go from the peak of the heating 
profile (see Fig. 4) at r ~ 12 km to R c , we find that r m is of 
order ~ 10 ms for L|J = 8.0 and r m ~ 1 second for L 5 ^ = 0.5. 
The steady-state approximation is only valid if r m , t s -C tj. For 
example, taking L v (t) oc f" 9 and high neutrino luminosities, M 
drops 10% in roughly tj. At these luminosities, is approx- 
imately 20-30 ms. Although both r m and t, are less than Td, 
they are not significantly so. We conclude that the steady-state 
assumption might be reasonably employed, but that caution is 
warranted. 



6. THE EVOLUTION OF PROTONEUTRON STAR WINDS 

With the eigenvalue problem solved and some of the system- 
atics in hand, in this section we explore possible evolutionary 
trajectories using the steady-state solutions. Beyond surveying 
the entire relevant parameter space, we endeavor to model the 
whole of the Kelvin-Helmholtz cooling phase, including radial 
contraction and the simultaneous evolution of the luminosity 
and average neutrino energy. 

Perhaps a second after core bounce, as the wind emerges, the 
protoneutron star atmosphere will be extended (R ~ 30-50 km) 
and perhaps highly luminous {L l ° l ~ 5 x 10 52 erg s" 1 ). As the 
neutron star cools it will contract quasi-hydrostatically. This 
may take as many as several seconds, depending upon the nu- 
clear equation of state. The average neutrino energies during 
contraction may increase, peak near the time at which R u set- 
tles, and then decrease roughly linearly in time (Pons et al. 
1999). The luminosity may decay quasi-exponentially or as a 
power law in time (Burrows and Lattimer 1986; Burrows 1988; 
Pons et al. 1999). 

Figure 10 shows the luminosity, radius, and average energy 
as a function of time for our evolutionary models. This pic- 
ture is merely schematic, but illustrates a representative sce- 
nario. In the following discussion we take L}° 1 oc f" 09 . A 
simple rescaling in time can be performed for other power- 
law indices or exponential luminosity decay. Two possible 
tracks for the time evolution of R u (t) are shown. The short 
dashed line is linear contraction such that R v (t = 0.4 s) = 20.3 
km and R u {t = Is) = 10 km. This is the evolutionary model, 
which we label as l R v {t) oc 1 —at'. For comparison, the dot- 
dashed line has R u (t) proportional to f" 1 / 3 . This model also 
has R u {t = 0.4s) = 20.3 km. At t ~ 3.2 seconds, R v = 10 km. 
The power 1/3 was obtained from a rough fit over approxi- 
mately one second of post-explosion evolution in a supernova 
model of S. Bruenn. In this calculation the supernova was sim- 
ulated in one-dimension artificially for a Woosley and Weaver 
(1995) 15 Mq progenitor starting from the collapse calcula- 
tions of Bruenn, De Nisco, & Mezzacappa (2001). We fo- 
cus on the model with R u {t) ocl—at instead of the model with 
R v {t) oc f" 1 / 3 because it reaches a more compact configuration 
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(i.e., maximum M/R v ) at earlier times, that is, with higher lu- 
minosity. As we explore in the next section, large M/R v , cou- 
pled with high luminosity and/or average neutrino energy gives 
short dynamical timescales and relatively high entropies, both 
potentially important for r-process nucleosynthesis in the pro- 
toneutron star context. 

In order to determine appropriate numbers for the ratios 

( £ V„ ) / ' ) / { £ fe ) > L Ve / L f„ > alld L Ve l L Ve WC SUn^d SU" 

pernova simulations (e.g., Mayle, Wilson, and Schramm 1987; 
Burrows, Hayes, and Fryxell 1995; Mezzacappa et al. 2001; 
Liebendorfer et al. 2001; Bruenn, De Nisco, and Mezzacappa 
200 1 ; Rampp and Janka 2000; S . Bruenn 200 1 , private commu- 
nication) and protoneutron star cooling calculations (Burrows 
and Lattimer 1986; Pons et al. 1999; J. Pons 2000, private com- 
munication). The common assumption of equipartition in lumi- 
nosity between the three neutrino species is generally not real- 
ized in these calculations. In fact, L v +L Vr +L„ t is usu- 
ally of order 50-60% of U° l . In addition, the ratio (e Pe ) /(e Ve ) 
ranges from 1.1 to 1.4 and L„jL Ve from 1.0 to 1.4. The ratio of 
(sv e ) to (e v ) also varies significantly. Like the fiducial mod- 
els presented in §5, we use Lp e to index our evolutionary mod- 
els in this section and set (s Vfl )/(e^ e ) = 1.6, (£p e )/(s Ve ) = 1-3, 
LvjL Vii = 1.4, and LyjL Ve = 1.3. At any time, all luminosities 
and average energies can then be computed from Fig. 10. We 
chose {Sv e )/(e^ e ) and L„jL Ve so as to accord with the litera- 
ture while also minimizing Y" (eq. 16), this being potentially 
favorable for r-process nucleosynthesis (but see §8.1). Note 
that with (£v e )/ (e Ve ) = 1 .3 and LyjL Ve = 1.3 and the magnitude 
of the average energies set by Fig. 10, Y° ~ 0.46 at early times. 
We feel these numbers are merely representative. We explore 
potentially important modifications to our prescription in §8.3. 

For each point along the evolutionary models represented in 
Fig. 10 we calculate the steady-state wind solution. We do this 
for a range of protoneutron star masses from 1.4-2.0 M Q . Neu- 
trino luminosity is always quoted as the local neutrino luminos- 
ity at the surface of the protoneutron star, not the luminosity at 
infinity. 

6.1. Results: Evolutionary Models 

Shown in Fig. 11 are evolutionary trajectories for M = 1.4, 
1.6, 1.8, and 2.0 M Q in the plane of r p versus asymptotic en- 
tropy, s a . Note that s a does not include contributions from (3- 
decays during nucleosynthesis. During r-process nucleosynthe- 
sis these processes may increase s a by > 10 units, depending on 
the dynamical timescale (Meyer and Brown 1997). However, 
because we post-process our wind models to obtain the nucle- 
osynthetic yield and include only alpha particles and free nucle- 
ons in the equation of state we employ in solving the eigenvalue 
problem (see §2), such an entropy increase is not included in 
s a . Figure 11 is analogous to Fig. 9, but for changing R v {t), 
(£„)(f), and L l ° l (t) using the evolution depicted in Fig. 10. The 
evolutionary trajectories labelled with 'R v oc I -at' are solid 
lines with small dots. For comparison, we also show the evo- 
lution for M = 1.4 M Q with R v (t) oc r 1 / 3 as a thin solid line 
without dots. The small dots on each evolutionary track are 
separate L„s. All tracks start with L 5 ^ = 8.0. The lowest lu- 
minosity shown on this plot is = 0.4 for each track. The 
time evolution for any mass begins with high luminosity, large 
R v , and, hence, low s a (<~ 50-70) and moderate r p (~ 9 ms). 
As R u gets smaller in the first second of evolution, the trajec- 
tories with R v oc 1 - at move to much higher s a and slightly 



smaller r p before they cease contraction at R v = 10 km. The 
s a reached at this luminosity is set by MjR v , with the 2.0 M Q 
model reaching s a ~ 150. Our evolutionary models contract 
from ~ 20 km to R v = 10 km in approximately 1 second. At 
this point in the evolution, R u is fixed and each track makes a 
sharp turn toward much longer r p and only moderately higher 
s a . This turnoff point is marked with a large open circle on 
each track and has L%\ = 3.4 and L l ° l = 1.57 x 10 52 erg s" 1 . At 
this point the trajectories join lines of constant R u , like those 
in Fig. 9. Due to the relatively slow contraction, the model 
with R u (t) oc f" 1 / 3 never exhibits such a sharp turn in the s a - t p 
plane and eventually joins the other 1.4 M Q evolutionary track 
at t p <~ 0.015 seconds, corresponding to t ~ 3.15 seconds and 
Lll = 1.3. Note that the turnoff point at L|J = 3.4 marks the 
point of minimum t p for each model with R v oc I -at. Table 2 
gives the global properties of our neutrino-driven wind models 

at Ll] = 8 - ' L U = 3 - 4 ' and L Ve = °- 4 ' mcmdm g the asymptotic 
electron fraction, Y". For comparison, we also include in Table 

2 the model with R u (t) cx r 1 / 3 at L 5 ^ = 3.4. 

Having solved for M at every point along these evolutionary 
tracks and assuming that L v (t) oc f" 09 , we calculate the total 
mass ejected in the wind as a function of time: 

M ej (f)= / M(t')dt'. (34) 
Jo 

Figure 12 shows this integral for all five of the models presented 
in Fig. 11. The dashed line shows M e j(f) for the model with 
R v (t) oc r 1 / 3 . As one would expect, because of the slower ra- 
dial contraction of the protoneutron star, this model ejects more 
matter than the corresponding trajectory with R v oc (I -at). In 
this case, the difference is about 30%. The small dots on each 
of the solid lines mark L v (t) for each model and correspond 
to the luminosity points on each track in Fig. 11. The large 
dots on each line mark L 5 ^ = 3.4, the luminosity at which each 
track in Fig. 11 reaches R v = 10 km and turns sharply. Note 
that the 2.0 M model ejects only ~6x 10" 5 M of mate- 
rial, about three times less than the 1.4 M Q model. Extrapo- 
lating the results of Fig. 12 we can compute the total M l ° l for 
t — ► oo. We can then compute, at any time, the mass yet to be 
ejected by the wind, AM e j(f) = M l ° l -M e j(f). Figure 13 shows 
AM e j(f) versus time for each track in Fig. 11. The lines and 
dots correspond with those in Fig. 12. In Fig. 1 1, we plot lines 
of constant log 10 [AM e j(f)] in units of M Q as dashed lines con- 
necting big dots on each of the four evolutionary trajectories 
with R u oc (1 -at). The thick dashed line on the far right side of 
the plot, labelled -6.0, is the line beyond which, only 10" 6 M Q 
will be ejected. 

In the calculations presented here, we have arbitrarily defined 
the point in time when the wind begins. The absolute magni- 
tude of M e j in Fig. 12 for each model is therefore also arbitrary. 
Only the ratios of these ejected masses or AM e j(f) for an indi- 
vidual trajectory are of real import. The -6.0 line in Fig. 1 1 
is of particular significance because if all (or most) supernovae 
produce r-process elements then the total yield per supernova 
must be 10" 5 - 10" 6 M (e.g., Qian 2000). Therefore, if an r- 
process epoch is to exist along any of the trajectories shown in 
Fig. 11, then it must begin at or before the line labelled -6.0 in 
order to eject sufficient mass. If r-processing begins to the right 
of this line, less than 10" 6 M Q will be ejected. For any L v (t) 
and R v {t), such a bound must exist. We have explored the po- 
sition of this boundary for a variety of relationships for L v (t). 
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Taking reasonable e-folding timescales (r) and L v {t) cx e~'l T 
the -6.0 line moves to even shorter r p . For slower power law 
decay, the boundary moves to longer r p . For M = 1.4 M Q and 
L v {t) cx t~ os it moves from r p ~ 0.07 to t p ~ 0.085 seconds. 
Although the wind may eventually evolve to arbitrarily long 
dynamical timescales, we conclude that the range of t p relevant 
for r-process nucleosynthesis is significantly constrained by the 
log 10 [AM e j] = -6 line in Fig. 11. In fact, this range is smaller 
than previous calculations suggest. In addition, note that even if 
a given wind model produces r-process elements, only a frac- 
tion of the total mass ejected during that r-processing epoch 
will be in r-process elements; much of the mass will remain in 
alpha particles. Conservatively, then, if transonic protoneutron 
star winds are the primary site for the r-process, this constraint 
on the amount of mass ejected per supernova implies that the 
epoch of r-process nucleosynthesis must occur for t p less than 
~ 0.085 seconds. 

For this range in t p , there is also only a relatively small 
range of s a available to the transonic protoneutron star wind. 
As evidenced by the calculations of Takahashi, Witti, & Janka 
(1994) and Qian & Woosley (1996) and borne out in Fig. 11, 
s a as large as ~ 400 is simply outside what can be obtained 
for reasonable dynamical timescales, even including the effects 
of general relativity. If we extrapolate the curves shown here 
to later times (lower luminosities), even the 2.0 M Q trajec- 
tory does not reach 400 until r p ~ 0.5 seconds. At this point 
M ~ 1.5 x 10~ 9 M© s" 1 . If r-processing could occur during 
these late stages, it would need to persist for many thousands 
of seconds to yield even 10~ 6 M of ejecta. Moreover, the 
survey calculations of Hoffman, Woosley, & Qian (1997) and 
Meyer & Brown (1997) show that with s a ~ 400 and r p ~ 0.5 
seconds, one requires Y" < 0.3 to achieve an appreciable third 
peak r-process. A Y° this low is extremely unlikely. As the 
protoneutron star cools, the neutrino luminosity and average 
neutrino energy will be correlated. Protoneutron star cooling 
calculations (Burrows and Lattimer 1986; Pons et al. 1999) in- 
dicate that the average neutrino energies will fall throughout 
these late evolutionary stages as the luminosity decreases. As 
the magnitude of (e Pe ) and (e Ve ) decrease, the asymptotic elec- 
tron fraction (eq. 16) must increase on account of the energy 
threshold for the reaction D e p — ► ne + (the 'threshold effect'). 
For the 1.4 M Q model, these constraints are even more severe. 
This track reaches s a <~ 400 only when r p is several seconds and 
M is of order 10" 11 M s" 1 . 

Note that for a given mass and t p , s a is 10-30 units higher in 
Fig. 9 than in Fig. 1 1 owing to the lower average neutrino ener- 
gies for a given luminosity in our fiducial models (§5) than in 
the evolutionary models we consider in this section. One might 
argue that by quickly decreasing (e u ) for all neutrino species 
with respect to the luminosity that the trajectory would move 
more quickly to higher s a , and, hence, be more likely to yield 
r-process ejecta. While such a change would certainly drive s a 
higher, it would also make r p increase faster, decrease M sig- 
nificantly, and the threshold effect would drive Y" higher. Thus, 
such a modification can only make the constraints tighter. 

We conclude that the late-time r-process as obtained in 
Woosley et al. (1994) is extremely unlikely in the context of a 
transonic wind. In essence, because s a is initially set by M/R v 
for a given model, once each trajectory reaches R u (t) = 10 km, 
the wind evolves quickly to much larger t p and only mod- 
estly higher asymptotic entropy. That is, for constant R u , s a 



is roughly proportional to r p 2 . By the time a transonic wind 
evolves to high entropy, the dynamical timescale is too long and 
Y" is too high to allow for a robust r-process. The slope of the 
trajectories in the s a - r p plane shown in Fig. 1 1 guarantee that 
if the wind enters a regime of very high entropy it does so with 
very large t p and minute M, so as to preclude any significant 
r-process yield. 

Instead, we propose an r-process epoch just a second or two 
after explosion, coinciding with the end of the protoneutron 
star contraction phase. In this scenario, the wind moves into 
an early-time r-processing regime in the s a - r p plane and then 
out of this regime at later times so that the constraint on M e j is 
satisfied. With this in mind, the behavior of the wind trajecto- 
ries during contraction, and particularly the point in each track 
where r p is at a minimum (L^ = 3.4), is suggestive and tantaliz- 
ing. As Hoffman, Woosley, & Qian (1997) noted, a small dy- 
namical timescale, even for only moderate entropies, can yield 
a successful r-process. More recently, Otsuki et al. (2000) have 
shown that a successful r-process can be realized in this con- 
text. For these reasons we turn our attention to an early-time, 
high luminosity, short- r p , and modest entropy r-process. For 
reference, in Fig. 1 1, we include long dashed lines of constant 
Y e , taken from the r-process survey calculations of Meyer & 
Brown (1997), above and to the left of which, for a given Y", 
production of the third r-process peak at A <~ 195 is assured. 
Caution is encouraged in taking these lines too seriously. They 
were computed along specific T and p trajectories for fixed Y° 
and dynamical timescale, r MB , given by eq. (33), which over- 
estimates drjdt by roughly a factor of two at these high lumi- 
nosities (see Fig. 8 and §5.1). These lines are only suggestive. 

7. NUCLEOSYNTHESIS IN THE EVOLUTIONARY 
PROTONEUTRON STAR WIND MODELS 

Any successful r-process site must do more than simply pro- 
duce nuclei with A > 195. Observations of r-nuclei in ultra 
metal -poor halo stars (notably, CS22892-052, HD 115444, and 
CS31082-001) show that the abundance pattern for A > 135 is 
nearly identical to the scaled solar r-process abundance (Sne- 
den et al. 1996; Burris et al. 2000; Westin et al. 2000; Hill et al. 
2001). Simply producing the platinum peak is no guarantee that 
the solar abundance distribution is reliably reproduced (Meyer 
and Brown 1997). In the wind scenario, particularly, one must 
construct the time-integrated yield as the neutrino luminosity 
decays and the global wind structure evolves. 

For each luminosity point on the 1 .4 M Q evolutionary trajec- 
tory with R v cx (1 - at) in Fig. 1 1 we obtained a unique M and 
velocity, temperature, and density profile. For each individual 
profile we compute the nucleosynthetic yield (Y) as a function 
of the atomic mass, A. With Y(A) and M at every point, assum- 
ing L„(f) cx f" 0,9 , we compute the weighted sum to get the total 
amount of ejected material at each A, M e j(A). For the 1.4 M Q 
model we find no significant nucleosynthesis beyond A ~ 100. 
In fact, most of the mass is concentrated at a peak in Sr, Y, and 
Zr. Inspecting the yield at each luminosity (time) reveals that 
when the protoneutron star has contracted to R v = 10 km, the 
point of minimum r p in Fig. 1 1 denoted by a large open cir- 
cle, the nucleosynthetic flow reaches a maximum in A. That 
is, all the points with t p > 0.006 seconds on the 1.4 M evo- 
lutionary track, even though they have higher entropy, produce 
lower average A ejecta. This can be understood simply: as the 
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dynamical timescale of the wind gets longer, more seed nuclei 
are formed. Hence, for a given Y" and s a , the neutron-to-seed 
ratio decreases (Hoffman, Woosley, and Qian 1996; Meyer and 
Brown 1997; Freiburghaus et al. 1999). Therefore, the point 
of minimum r p , when MjR v reaches a maximum, affords the 
best possibility for a robust r-process. For all protoneutron star 
masses, the evolutionary models with R v oc (1 -at) turn sharply 
at L|' =3.4 and for t p between 6 and 7 ms. Unfortunately, al- 
though the Lp' = 3 .4 point produced the highest average A ejecta 
of any other luminosity along the 1.4 M Q track in Fig. 1 1, the 
nucleosynthesis did not even proceed to the second abundance 
peak. 

We also calculated the nucleosynthetic yield for the 1.6, 1.8, 
and 2.0 M Q trajectories at Lp 1 = 3.4, assuming that these points 
of minimum r p would also yield the highest average A ejecta 
of any of the points in a given mass trajectory. None success- 
fully generated nucleosynthesis beyond the second r-process 
abundance peak. Even the 2.0 M Q model, which has s a ~ 151 
and t p ~ 0.0068 seconds for L 5 J = 3.4, did not proceed beyond 
A ~ 135. Hence, for the Y"s derived and the time evolution 
we have assumed, we fail to produce viable r-process nucle- 
osynthesis in any of our evolutionary transonic protoneutron 
star wind models. 

8. WHAT IS TO BE DONE? 

We have already ruled out the possibility of a late-time r- 
processing epoch at long dynamical timescales and high en- 
tropies (s a > 400) using constraints on the amount of material 
ejected, the slope of the evolutionary tracks in the s a - t p plane, 
and the inexorable rise in Y" as the protoneutron star cools. 

We are left wondering what reasonable modifications might 
generically yield third-peak r-process nucleosynthesis for a 
canonical protoneutron star with 1.4 M Q and R v = 10 km. In 
this section we further explore the viability of the early-time r- 
process. We consider the likely range of Y", s a , and r p accessi- 
ble to transonic neutrino-driven winds, and present the physical 
conditions we require for production of the both the second and 
third r-process abundance peaks. 

8.1. The Asymptotic Electron Fraction: Y" 

One might argue that Y° is simply too high in these winds 
to yield successful nucleosynthesis. There are several im- 
portant points in this regard. First, our evolutionary models, 
which failed to produce nucleosynthesis beyond A ~ 135, all 
had Y" — 0.46 at L 5 ^ = 3.4. This Y° favors the formation of 
A <~ 90 nuclei and produces many seed nuclei, thus decreas- 
ing the neutron-to-seed ratio for a given entropy and dynamical 
timescale. Second, both Woosley et al. (1994) and Wanajo et 
al. (2000) obtained unacceptably large over-productions of nu- 
clei near A ~ 90 in the early phase of their wind calculations, at 
high luminosity and low entropy. Third, Hoffman et al. (1996) 
find that this overproduction problem is solved if Y" > 0.485. 
These three points together imply that if the r-process occurs 
generically in protoneutron star winds then Y° must be either 
less than 0.40 or greater than 0.48 to avoid the overproduction 
problem at A ~ 90. Naively, it might seem that Y" < 0.40 is 
favored because this would naturally increase the neutron-to- 
seed ratio by increasing the number fraction of neutrons. How- 
ever, there are several reasons why Y" > 0.48 might actually 



be viable. First, for Y e a > 0.485 Hoffman et al. (1996) found 
that some interesting p-process elements were produced, which 
were previously unaccounted for (e.g., 92 Mo). Second, the most 
detailed transport studies done to date (Mezzacappa et al. 2001; 
Liebendorfer et al. 2001; Rampp and Janka 2000) indicate that 
(Sv e )/ (£v e ) — 1.1 — 1.2 and L v jL Ve ~ 1.1. Depending on the 
magnitude of (e Pf ) and (e Ve ), these numbers put Y° > 0.48, as 
per eq. (16). Third, Hoffman, Woosley, & Qian (1997) find that 
as Y" increases from 0.48 to <~ 0.495 the requisite entropy for 
third-peak production actually decreases for fixed dynamical 
timescale. The last point implies that having a high Y" might 
slightly relieve the constraints on s a set by an early-time, high- 
luminosity r-process. The fact that some /^-process elements 
might also be produced in a high Y" environment is attractive. 
Together, we feel that the above points make it plausible that 
Y" > 0.48 in protoneutron star winds. Such a conclusion, con- 
strains the three-dimensional space s a — T p - Y" significantly. 

Of course, having Y" ~ 0.30-0.35 might also cure the over- 
production problem at A ~ 90 while increasing the neutron-to- 
seed ratio dramatically, so as to allow for third-peak produc- 
tion at the entropies and timescales obtained for the 1.4 M Q 
evolutionary model in Fig. 11. However, to attain Y" < 0.35 
one requires UjL Ve > 1.55 for (e Se )/<£„,) = 20 MeV/12 MeV. 
Such conditions would be extreme in light of the detailed trans- 
port calculations carried out to date. However, if Y" does reach 
values this low in profiles like those in Fig. 11, there are two 
constraints worth pointing out. The first is that Y" must evolve 
with the neutrino luminosities and average energies so that very 
little mass is left to be ejected by the time Y" increases to 0.40. 
Otherwise, the same overproduction problems at A ~ 90 may 
occur. The second constraint is that if Y" is sufficiently low 
to guarantee third-peak production, it must eject not more than 
<~ 10~ 5 M Q of r-process material per supernova. 

We conclude that Y° may be >0.47 in protoneutron star 
winds. This follows from the fact that Y" below 0.40 is very 
unlikely and if Y° is in the range 0.40-0.46, models suffer from 
overproduction of A <~ 90 nuclei. With this in mind, in the next 
section we consider modifications to our transonic wind mod- 
els that might increase the entropy or decrease the dynamical 
timescale. 

8.2. Entropy and Dynamical Timescale 

One might choose to increase s a by changing the bulk pro- 
toneutron star characteristics. Increasing the ratio MjR v in- 
creases s a significantly with only modest decreases in t p . How- 
ever, this ratio cannot be increased arbitrarily. M > 1 .5 M Q 
may be disfavored in light of neutron star binary observations 
(Arzoumanian 1995) and R u < 9 km seems unlikely due to con- 
straints on the high density nuclear equation of state (e.g., Lat- 
timer and Prakash 2001). In order to explore this, however, 
we varied M and R u at L 5 ^ = 3.4 in our evolutionary models 
with R v oc 1 —at. In Table 3, we summarize these results for 
M =1.8 M , 1.6 M Q , and 1.4 M . These models should be 
compared with the models with L|J = 3.4 in Table 2. Unfortu- 
nately, for Y" ~ 0.46 for each model, we did not obtain third- 
peak r-process nucleosynthesis. Increasing Y° artificially in our 
nucleosynthesis calculations to 0.48, the neutron-to-seed ratio 
stays low and we fail to generate r-process elements beyond the 
second peak. 

Although reasonable increases in M/R v are favorable for the 
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r-process, they are insufficient for strong third-peak nucleosyn- 
thesis. Therefore, in an effort to obtain a robust r-process in a 
canonical 1.4 M Q , R v =10 km protoneutron star, we are only 
left with the option of modifications to the energy deposition 
profile. 

8.3. Possible Alterations to q 

Qian & Woosley (1996) showed that an artificial energy 
source at radii between 20 km and 30 km, beyond the peak in 
q could substantially increase the entropy and decrease the dy- 
namical timescale. In fact, any extra energy source that broad- 
ens the energy deposition profile, thus increasing q in a region 
of low mass density, increases s a and decreases t p . 

We noted in §6 the difference in entropy, for a given t p , be- 
tween the fiducial tracks in Fig. 9 and the evolutionary tracks 
in Fig. 11. Comparing the M =1.6 M Q tracks on both plots, 
the difference in s a between the two at r p = 0.02 seconds is 
<~12 units. The fiducial model (with higher s a ) has a total neu- 
trino luminosity almost twice that of the evolutionary model, 
but its average neutrino energies are more than 35% lower. 
The increase in entropy is caused by an interplay between 
the charged-current heating rate (q cc ) and the neutrino-electron 
scattering heating rate (<7„ e ). The former is proportional to 
L v e {£ 2 v e ) +L Ve {e1 t ), while q Vie is proportional to J2i L ^( £ ^,)- 
Clearly, for fixed neutrino luminosities, as the average neutrino 
energy drops, q Vie becomes more important relative to q cc . Be- 
cause q cc — > as the alpha fraction increases (see Figs. 6 and 
7), the fact that q v<e increases in importance effectively broad- 
ens the energy deposition profile, thus increasing the entropy. 
Although the heating rate due to v(Di annihilation peaks close 
to the protoneutron star surface, it also contributes to the to- 
tal energy deposition rate at radii larger than where q cc — > 0. 
Because q v . . oc X^-^v( £ ^)> f° r fixed neutrino luminosity and 
decreasing average neutrino energy, this process also becomes 
more important relative to q cc , thereby enhancing the effect on 
s a . Although the total effect here is relatively small, at short 
t p any increases in s a are of potential significance. Small av- 
erage energies coupled with higher luminosities are one way to 
achieve moderately higher s a and shorter t p . Note that due to 
the threshold effect (eq. 16) in the charged-current reactions, if 
one decreases (e Pf ) and (e Ve ), Y" will increase, and any poten- 
tial gains in s a might be mitigated. However, as we discussed in 
§8.1, Fig. 10 of Hoffman, Woosley, & Qian (1997) shows that 
the entropy required for third-peak nucleosynthesis actually de- 
creases for high Y". 

Similarly, one might also increase (e„ ) and L u ^ relative to 
the same quantities for the electron and anti-electron types. 
Since the v^- and j/ T -type neutrinos do not participate in the 
charged-current reactions, any increase in their luminosity or 
average energy effectively increases the importance of q Vi( , and 
q Vi Vi with respect to q cc . As we noted in §6, the evolution 
of luminosity and energy shown in Fig. 10 is only sugges- 
tive. For this reason we explored modifications to the ra- 
tios (£„„)/<£*.) = 1.6, (ep, )/<£„.) = 1.3, L P JL^ = 1.4, and 
Ljy e I L Ve = 1.3. For our extreme 2.0 M Q evolutionary model with 
L 5 J =3.4 (see Table 2 for comparison), we set L Ve = L Pe = L v 
and (e Ve ) = 15 MeV, {e Pe } = 22 MeV, and (e V)l ) = 34 MeV. This 
increased Y" from 0.469 to 0.484, increased Q by more than a 
factor of three, decreased s a by 17 units to ~134, and decreased 
T p by 40% to 4. 1 ms. These modifications were insufficient to 
produce third-peak nucleosynthesis. 



Following Qian & Woosley (1996), we artificially increased 
q(r) in the region 20 km < r < 50 km for our 1 .4 M Q evolution- 
ary model at L%\ = 3.4 so that Q (eq. 14) went from 1 .21 x 10 49 
erg s" 1 to 1.33 x 10 49 erg s" 1 , an increase of 10%. Because 
Q is the volume integral of pq and p is small at these radii 
(10 6 — 10 7 g cm -3 ), q must be enhanced substantially in order 
to affect a 10% change in Q. With this extra energy deposition 
we found that maximum increases in s a and decreases in r p , de- 
pending on the degree of augmentation of q as a function of r, 
were 17 units and 50%, respectively. M increased by just 8%. 
We made the same sort of modification to the 2.0 M Q model 
with L Vt = L„ e = L Vfi and (e Ve ) = 15 MeV, (e Pe ) = 22 MeV, and 
{e v ) = 34 MeV. In this case, we increased Q by 6% and found 
that s a was increased from 134 to 150 and that t p decreased 
from 4. 1 ms to 2.6 ms. That t p can decrease so significantly as 
a result of < 10% changes in Q demonstrates the importance of 
conducting a full transport study in this context. We save such 
an investigation for a future work, but emphasize that the shape 
of the energy deposition profile may be the final arbiter in de- 
termining the true potential of this site as the seat of r-process 
nucleosynthesis. 

8.4. The Early-Time r-Process 

Although we have described the general physics of protoneu- 
tron star winds, the resulting nucleosynthesis, and modifica- 
tions to our models that might enhance the wind's entropy 
and decrease its dynamical timescale, none of the models we 
have presented so far produces a robust, third-peak r-process. 
However, setting M = 2.0 M Q , R v = 9 km, L Pc = L Ve = L„ M = 
8.0 x 10 51 erg s~\ and (e Ve ) = 14.5 MeV, (e p J = 22 MeV, 
(e Vtl ) = 34 MeV, and adding an artificial heating source be- 
tween 20 km and 50 km that increases Q by 4%, we derive a 
wind with s a ~ 150, t p ~ 1.3 ms, and Y" ~ 0.477. This ex- 
tremely compact and luminous protoneutron star yields a wind 
profile that produces third-peak r-process nucleosynthesis. We 
have found that for s a ~ 150 and 0.47 < Y e a < 0.495, we re- 
quire t p < 1.3 ms in order to generate a significant A <~ 195 
yield. These are the necessary conditions we derive from our 
general-relativistic wind models for an early-time protoneutron 
star r-process epoch. If a wind trajectory, like those in Fig. 1 1 
were to pass into the region s a <~ 150 and t p < 1 .3 ms with Y° 
less than 0.50, some third-peak material would be produced. 
Note that Otsuki et al. (2000) attained third-peak nucleosyn- 
thesis for similar wind conditions. Artificially setting Y° equal 
to 0.40, their model with s a ~ 140 and t p ~ 1 .2 ms successfully 
produced abundance peaks at A <~ 130 and A ~ 195. 

Although a full nucleosynthesis survey, using real wind pro- 
files, is required to map the space 100 < s a < 200, t p < 1 .5 ms, 
and 0.46 < Y e < 0.50, we note some features of potential im- 
portance for the short r p , early-time r-process. First, as we de- 
crease Y" from 0.495 to 0.47, for a given t p and s a , the neutron- 
to-seed ratio stays roughly constant. Hence, the ratio of the 
abundance yield at A <~ 130 to that at A <~ 195 is relatively in- 
sensitive to Yg. Second, r-process nucleosynthesis at very short 
timescales and high electron fractions is possible because the 
number of seed nuclei formed is very small. As a consequence, 
we expect the nucleosynthetic yield in this regime to be sen- 
sitive to changes in the input nuclear physics and, in particu- 
lar, the three-body reactions important in seed nuclei formation 
(e.g., 4 He(a«,7) 9 Be; Kajino et al. 2001). 



15 



Our requisite conditions for third-peak nucleosynthesis, t p < 
1.3 ms, s a ~ 150, and high Y", disfavor r-process nucleosynthe- 
sis generically in neutrino-driven winds from neutron stars with 
M = 1 .4 M Q . Our results in Table 3 indicate that even for R u = 8 
km, s a is 35 units too small and r p is a factor of about three too 
long for the r-process to proceed to the third abundance peak. 
We find these gaps in entropy and timescale very difficult to 
bridge. For M = 1 .4 M Q we require R v < 6.5 km to reach this 
s a and t p . It is unlikely that any high-density nuclear equation 
of state could accommodate such a small radius. Even taking 
M = 1 .6 M Q and R v =9 km, without invoking an artificial heat- 
ing source, we fail to reach t p < 1 .3 ms and s a ~ 150. Although 
the importance of the distribution of energy deposition and ex- 
tra heating sources should be borne in mind, our unmodified 
wind models require a very massive and highly luminous pro- 
toneutron star with small radius. Indeed, considering the fact 
that our successful wind models originate from neutron stars 
with M > 2.0 M and R v < 9 km, we are forced to consider 
the possibility that the primary site for the r-process is not a 
protoneutron star at all. A neutrino-driven outflow generated 
near the event horizon of a black hole might bear many of the 
characteristics of our successful protoneutron star models. Per- 
haps very short timescale outflows or jets originating from the 
compact inner accretion disk created in the collapsar models of 
MacFadyen & Woosley (1999) attain the necessary conditions 
for r-process nucleosynthesis. Such outflow models would ben- 
efit by being generated in a region with high M/R v , without 
being subject to the constraints imposed on neutron stars by the 
nuclear equation of state. 

9. DISCUSSION 

We have constructed a robust and efficient algorithm for solv- 
ing the neutrino-driven protoneutron star wind problem using 
both general relativity and Newtonian gravity. We employed 
physical boundary conditions for the transonic wind, a well- 
motivated neutrino energy deposition function, and an equation 
of state suited to this problem. For the first time, we included 
the differential equation for the evolution of Y e in radius and 
the proper sonic point boundary condition. Using this computa- 
tional tool, we studied the structure and systematics of neutrino- 
driven winds with an eye toward assessing the suitability of this 
site for r-process nucleosynthesis. We have examined a wide 
range of protoneutron star radii, masses, and neutrino spectral 
characteristics. By positing an expression for L l ° l (t), we have 
modeled potential contraction and cooling scenarios that might 
exist in Nature and calculated the total mass ejected for the cor- 
responding evolutionary trajectories. Employing general rel- 
ativistic hydrodynamics, we find significant enhancements in 
the asymptotic entropy of the wind, in good agreement with the 
post-Newtonian models of Qian & Woosley (1996), the ana- 
lytic approximations of Cardall & Fuller (1997), and the work 
of Otsuki et al. (2000). In addition, we find that modest modifi- 
cations to the net energy deposition rate can markedly improve 
the conditions for r-process nucleosynthesis. Indeed, we feel 
that changes in the profile of energy deposition represent the 
most viable alterations to our models, which might lead to ro- 
bust r-process nucleosynthesis in the protoneutron star context. 

Our results indicate that only an early-time epoch of r- 
process nucleosynthesis at high L v , small r p , and modest en- 
tropy is possible. A late-time r-process, at very high entropy 
(> 300), long t p , and low L v is not viable. There are several 



components to this argument. As the luminosity of the pro- 
toneutron star decays, both the asymptotic entropy and dynam- 
ical timescale of the wind increase. The former is conducive to 
r-process nucleosynthesis. The latter is not. Hence, the inex- 
orable rise in both compete. Fundamentally, for the transonic 
wind, we find that the asymptotic entropy does not increase fast 
enough to compensate for the deleterious rise in r p . In addition, 
as the luminosity decays and r p increases, M decreases. For ex- 
ample, our 1 .4 M Q evolutionary model only reaches s a <~ 300 
when M~5x 10" 10 M Q s" 1 and t p ~ 5 seconds. Clearly, 
even if the r-process could exist under these conditions, such an 
epoch would have to continue for many thousands of seconds 
to produce even 10" 6 M of r-process ejecta. Finally, the con- 
tinued rise in Y" at late times as (s Ve ) and (e Pe ) decrease (due to 
the energy threshold for the reaction v e p — ► ne + , the threshold 
effect) also argues against a high-s fl , long-r p r-process. 

For these reasons, we conclude that if r-process nucleosyn- 
thesis occurs in protoneutron star winds, it must occur at early 
times, at or just after the moment when R u reaches a minimum. 
Our 1 .4 M Q evolutionary trajectory with R u oc \-atm Fig. 1 1 
does not attain sufficiently high entropies and short timescales 
for successful third peak r-process nucleosynthesis. We have 
calculated the nucleosynthetic yield at every luminosity (time) 
in this trajectory and at no point does the resulting nucleosyn- 
thesis go beyond A <~ 100. Interestingly, however, we find that 
the luminosity point that yields the highest average A ejecta cor- 
responds to the point in time where R u reaches a minimum. 
This point is also a minimum in r p for the trajectory. We calcu- 
lated the nucleosynthesis at this same luminosity point for each 
mass in Fig. 11. Even for the 2.0 M model with s a ~ 150, we 
did not obtain nucleosynthesis beyond the second abundance 
peak. As evidenced by the survey calculations of Hoffman, 
Woosley, & Qian (1997) and Meyer & Brown (1997) and the 
wind calculations of Otsuki et al. (2000), these models are out- 
side a regime of successful third-peak nucleosynthesis. How- 
ever, in §8.3 we have shown that reasonable modifications to 
the spectral character of the neutrinos and the energy deposi- 
tion function might conceivably shorten t p sufficiently for the 
r-process to proceed in some of these models. In §8.4, we found 
that winds with s a ~ 150, t p < 1.3 ms, and 0.47 < Y" < 0.495 
can generate third-peak r-process elements. 

If transonic protoneutron star winds are the primary site for 
r-process nucleosynthesis, then a successful r-process wind 
model must enter this s a —T p —Y" regime. The wind, starting just 
after re-ignition of the supernova shock, begins with large R v 
and L'° l and, hence, low entropy (s a <~ 50) and short dynamical 
timescales (r p is several ms). In order to avoid overproduction 
of A <~ 90 nuclei, Y" is high (> 0.48) during this low entropy 
contraction phase (see §8.1; Hoffman et al. 1996). As the pro- 
toneutron star contracts, it moves to much higher entropy and 
shorter t p . Just as R„ reaches a minimum and the protoneutron 
star is at its most compact, s a is sufficiently high (<~ 150), and 
Tp is sufficiently short (< 1.3 ms) to guarantee successful third 
peak r-process nucleosynthesis. This epoch does not persist. 
Because R v is now constant in time the wind evolves quickly 
along trajectories like those in Figs. 9 and 1 1 with s a oc r p 2 to 
much longer timescales and only moderately larger s a . This, 
coupled with the rise in Y° due to the threshold effect, effec- 
tively shuts off the r-process so that no more than <~ 10" 5 M Q 
is ejected. We emphasize that this scenario requires the wind to 
move into and then out of an r-processing regime in the space 
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of s a -r p -Y e a . 

Such a picture is provocative, but not yet convincing. Simply 
obtaining a wind solution that has proper s a , r p , and Y" to guar- 
antee production of the third peak is hardly sufficient to explain 
the remarkable agreement between the r-element abundances 
with atomic masses at and beyond barium in ultra-metal-poor 
halo stars and the observed solar r-process inventory. It is dif- 
ficult to understand how such a scenario might consistently re- 
produce the barium abundance, all the lanthanides, the platinum 
peak, and the actinides (Cayrel et al. 2001; Hill et al. 2001). In 
addition, while the data for these halo stars show remarkable 
consistency with the solar abundances above A <~ 135, below 
this mass they are markedly inconsistent and there is signifi- 
cant star-to-star scatter. Perhaps some subset of all supernovae 
account for the region below A <~ 135 and never undergo a vig- 
orous r-process. Perhaps others do obtain the required s a , t p , 
and Y" and account for the full range of nuclides, including 
the first abundance peak and proceeding to uranium (Wasser- 
burg and Qian 2000; Qian and Wasserburg 2000; Sneden et al. 
1996; Burris et al. 2000). As we demonstrate in §6 with our 
R v oc f" 1 / 3 model, slow radial contraction of the protoneutron 
star may preclude any significant r-process yield, as the wind 
would then never enter a regime of short dynamical timescale 
with s a ~ 150. 

The supernova progenitor structure might be important in 
this regard. The two-dimensional calculations of Burrows, 
Hayes, & Fryxell (1995) and Janka & Muller (1995) indicate 
that a transonic protoneutron star wind can form just tenths of 
seconds after the successful re-ignition of the supernova shock. 
The pressure of the wind is sufficient for it to emerge into the 
expanding supernova ejecta. However, Janka & Muller (1995) 
found that as the supernova shock passes through the Si-O in- 
terface in their one-dimensional 15 M Q progenitor it causes 
a strong reverse shock that slows the wind expansion from a 
v~2x 10 9 cm s" 1 to a few times 10 8 cm s" 1 . It is possible 
that a termination or reverse shock might generally disrupt the 
transonic wind as it propagates toward the protoneutron star. 
Exactly how far in radius the reverse shock propagates will be 
a function of the hydrodynamical power of both the wind and 
the reverse shock as the neutrino luminosity decays, each be- 
ing functions of the progenitor structure. With sufficient power, 
the reverse shock may continue to the sonic point. This would 
put the whole region between the protoneutron star surface and 
the supernova shock in sonic contact, thereby converting a tran- 
sonic wind into a subsonic breeze. We have conducted prelim- 
inary hydrodynamical calculations, which suggest this could 
occur in certain circumstances. Steady-state solutions to the 
wind problem can also be formulated in this context, but with 
an outer boundary pressure set by conditions at the supernova 
shock. Qian & Woosley (1996) explored the effects of an ex- 
ternal boundary pressure on their wind models. This increased 
the wind entropy by just 11 units, but increased the dynamical 
timescale by more than 60%. Such a change would be detri- 
mental to an early-time, short- r p r-process. However, it may be 
that the reverse shock does not have sufficient power to disrupt 
the wind interior to the sonic point. In the steady-state, across 
the shock boundary, the velocity will decrease, the density will 
increase so as to maintain M, and both the temperature and en- 
tropy will increase. 

In order to test the effects of a termination or reverse shock 
on the nucleosynthesis, we inserted a shock by hand at a radius 



(r sh ) of 4000 km in our wind model with M = 2.0 M and = 9 
km (see §8.4), far outside the sonic point (r ~ 180 km). At r s h 
the matter velocity was 5.1 x 10 9 cm s" 1 , p was 17 g cm" 3 , and 
T was approximately 0.023 MeV. Using the Rankine-Hugoniot 
shockjump conditions, we estimate that v ~ 7v', p ~ p' /7, and 
T <~ T'/2, where unprimed quantities are for the wind just be- 
fore the shock (r < r s h), and primed quantities are for the flow 
just after the shock (r > r s h). These conditions increase s a by 
about 10 units in the post-shock region and increase r p signifi- 
cantly, due to the sudden decrease in v. This had subtle, but po- 
tentially significant effects on the resulting r-process yield. To 
appreciate this, one must understand that without the slowing of 
the fluid trajectory by passage through the reverse shock, the r- 
process freeze-out in our wind models occurs for temperatures 
below 0.01 MeV, because only then are beta decays along the 
r-process path fast enough to compete with the rapid material 
expansion. By contrast, if the material slows (and reheats) by 
passage through the shock, the r-process freeze-out happens at 
higher temperatures, typically near 0.05 MeV. Although the av- 
erage number of neutrons captured per seed nucleus is the same 
for the shocked and unshocked (but otherwise identical) trajec- 
tories, the distribution of those neutron captures is different. In 
particular, for the case considered here, the shocked trajectory 
had a factor of three larger yield at A~495. The reason is that, 
when the trajectory slowed and reheated by shock passage, the 
nuclear flow changed and allowed more nuclei to leak out of 
the N=82 closed shell and proceed up to N=126 (A^195) at the 
expense of flow from the N=50 closed shell (A^80) to N=130. 
An additional interesting effect was that, unlike the unshocked 
trajectory, the shocked trajectory showed evidence of formation 
of a rare-earth element peak at A^165. Surman et al. (1997) ar- 
gued that this peak forms during freeze-out as the r-process path 
rapidly moves through the Z^60, N~404 region of somewhat 
enhanced nuclear stability in the nuclide chart. In the winds 
we study here, the shocked trajectories favor such a freeze-out 
while the unshocked trajectories do not. We conclude that the 
finer details of the r-process abundance curve may depend in 
interesting ways on the location and strength of a termination 
or reverse shock. 

These hydrodynamical issues are part of the larger ques- 
tion of fallback in Type-II supernovae. It is possible that the 
most massive supernova progenitors, with their extended hy- 
drogen envelopes and dense core structures, experience signifi- 
cant fallback over timescales much longer than the cooling time 
(Chevalier 1989; Woosley and Weaver 1995). Even if the wind 
were able to emerge from the neutron star for 10-20 seconds 
after explosion, it might not have sufficient power to overcome 
extended fallback over minutes and hours. Even without a large 
overlying hydrogen envelope (Type lb, Ic), the neutrino-driven 
wind may be hindered by any progenitor with a large inner core 
and outer core binding energy. Therefore, we speculate that 
an early-time r-process, unencumbered by fallback or reverse 
shocks, is most likely in less massive Type-II, -lb, or -Ic super- 
nova progenitors. Accretion-induced collapse may offer even 
more potential in this context for a fully developed, early-time 
transonic wind as, in this case, there is no overlying mantle to 
impede the wind's emergence (Fryer et al. 1999). 

10. CONCLUSIONS 

Our major conclusions are the following: 
• For a given mass outflow rate, we find a significantly 
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shorter dynamical timescale (t p ) than indicated by 
many previous investigations. 

• Because the temperature, density, and velocity 
gradients often used in defining a dynamical timescale 
evolve on a mass element as it moves away from the 
protoneutron star, employing a unique dynamical 
timescale to characterize the wind is not recommended. 

• For a given protoneutron star radius, the asymptotic 
entropy (s a ) is proportional to r a p 2 . 

• A late-time, high entropy (s a > 300), long timescale 
neutrino-driven wind, is not a viable astrophysical site 
for r-process nucleosynthesis. Although the wind may 
eventually realize very high entropy, the mass outflow 
rate will be too small and the dynamical timescale 
too long at such an epoch to account for the galactic 
r-process abundance. 

• The third-peak r-process elements can be produced in 
significant abundance in protoneutron star winds only 
at early times, at modest entropy (s a <~ 150), and very 
short dynamical timescale (r p <~ 1 ms). 

• Winds originating from protoneutron stars with mass 

1 .4 M Q and radius R v = 10 km do not produce elements 
beyond A <~ 100 at any time during wind evolution. 

• We derive that third-peak r-process arises naturally in 
the context of spherical, transonic protoneutron star 
winds only in the unlikely case of protoneutron stars 
with M > 2.0 M and R v < 9 km. 

• Shocks in the protoneutron star wind, exterior to 
the sonic point, caused by the wind's interaction 
with the inner supernova can significantly influence 
the third-peak and rare-earth r-process element 
abundances. 

All of these conclusions hold for the generally high Y/s 
(0.46-0.49) derived in this work. Only in the unlikely case 
that Y e in the wind is < 0.35 are our conclusions dramati- 
cally altered. In this case, third-peak r-process nucleosynthesis 
might be obtained in our transonic models for protoneutron star 
masses less than 2.0 M Q and radii greater than 9 km, at early 
times in the wind evolution. The uncertainty in the spectral 
characteristics of the electron and anti-electron neutrinos in de- 
termining the asymptotic electron fraction during wind forma- 
tion and evolution is primary on our list of unresolved issues. 
Others include: 

• the character of the energy deposition profile as 
obtained from detailed neutrino transport in a 
self-consistent calculation, 

• the high-density nuclear equation of state, which would 
elucidate the range of protoneutron star masses and 
radii relevant, 



• the hydrodynamical interaction of the wind and shocks 
in the expanding post-supernova environment, 

• the effects of rotation and magnetic fields, 

• exactly how, even given successful third-peak 
nucleosynthesis for a given model, the universality of 
the observed r-process distribution for A > 135 can be 
accounted for generically by the progenitor-dependent 
parameter space of neutrino-driven winds. 

These last points leave the prospect of the r-process in this 
context an open question. We conclude from this study that if 
the r-process occurs in protoneutron star winds, it most likely 
occurs at early times after the preceding supernova, in winds 
with very short dynamical timescales (< 1.3 ms), moderate en- 
tropies (<~ 150), and, possibly, high electron fractions (0.47< 
Y" < 0.495). In contrast to the early-time scenario, because 
protoneutron star winds only enter a high entropy (s a >300) 
regime with very low M (< 10~ 9 M Q s" 1 ) and extremely long 
timescales (r p > seconds), a late-time r-process is simply not 
viable. Conditions suitable for an early-time r-process are re- 
alized in our models only by very compact protoneutron stars 
with M = 2.0 M Q and R u < 9 km. The conditions neces- 
sary for third-peak r-process nucleosynthesis are not realized 
in neutrino-driven transonic winds from canonical neutron stars 
with M = 1 .4 M Q and R v = 10 km. In fact, for this neutron star 
mass, we require R v < 6.5 km. Although such neutron star 
masses and radii are not entirely excluded by current high den- 
sity equations of state, such radii seem unlikely to obtain in the 
early post-supernova phase. The short-timescale jet outflows 
from the dense inner accretion disk around a black hole formed 
in the collapsar models of MacFadyen & Woosley ( 1 999) might 
attain the necessary entropies and timescales for the r-process, 
since in that context MjR v can be significantly larger than in 
the protoneutron star context. Importantly, it should be noted 
that we consistently produce r-process nucleosynthesis below 
A ~ 135. Perhaps protoneutron stars of canonical mass and 
radius (1.4 M Q , 10 km) produce elements in this mass range 
generically, thus accounting (due to progenitor structure and 
temporal characteristics of the neutrino spectrum) for the vari- 
ations in abundance observed in these elements in ultra-metal- 
poor halo stars. 
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Table 1 

Fiducial Wind Models: 1.4M f 



L % 


L'° l (10 51 erg s ') 


M (M s ') 


Q (10 48 erg s ') 


fmcch (10 48 erg s ') 


t p (ms) 




8.0 


37.0 


9.05 x ur 5 


35.1 


0.848 


3.68 


83.9 


DC 

b 


37.0 


2.70 x ur 4 


75.1 


2.54 


3.28 


68.2 


7.0 


32.4 


6.56 x Hr 5 


25.1 


0.561 


4.24 


86.4 


6.0 


27.6 


4.47 x Hr 5 


17.1 


0.346 


5.01 


89.5 


5.0 


23.1 


2.84 x Hr 5 


10.9 


0.194 


6.14 


93.4 


4.0 


18.5 


i.63 x ur 5 


6.28 


0.0929 


7.93 


98.4 


3.0 


13.9 


8.03 x ur" 


3.09 


0.0346 


11.29 


105.4 


2.0 


9.25 


3.74 x ur" 


1.14 


8.59 x Hr 3 


19.19 


116.2 


1.0 


4.63 


5.44 x 10~ 7 


0.211 


8.24 x Hr 4 


49.54 


137.6 


0.6 


2.78 


i.58 x io~ 7 


0.062 


1.50 x 10" 4 


100.9 


156.3 



Newtonian calculation 



Table 2 

Evolutionary Wind Models: Ll\ = 8.0, 3.4, and 0.4 



Mass (Mq) 


L- 1 


L™(10 51 ergs" 1 ) 


R u (km) 


M(M s- 1 ) 


GdO^ergs- 1 ) 


fmcchUO 48 ergs" 1 ) 


t p (ms) 




Y° 


2.0 
1.8 
1.6 
1.4 


8.0 
8.0 
8.0 
8.0 


37.0 
37.0 
37.0 
37.0 


20.3 
20.3 
20.3 
20.3 


2.19 x ur 4 

3.03 x Hr 4 
4.38 x Hr 4 
6.69 x Hr 4 


53.3 
67.2 
83.8 
108.2 


1.27 
1.67 

2.27 
3.25 


10.45 
9.87 
9.27 
8.63 


67.08 
59.05 
51.39 
44.02 


0.467 
0.465 
0.463 
0.460 


2.0 
1.8 
1.6 
1.4 
1.4° 


3.4 
3.4 
3.4 
3.4 
3.4 


15.7 
15.7 
15.7 
15.7 
15.7 


10.0 
10.0 
10.0 
10.0 
14.7 


9.80 x ur 6 
i.39 x ur 5 
2.04 x ur 5 
3.14 x ur 5 

9.03 x Hr 4 


5.94 
7.34 
9.29 
12.1 
21.8 


0.0826 
0.108 
0.148 
0.213 
0.424 


6.83 
6.78 
6.56 
6.21 
9.98 


151.61 
129.43 
109.70 
91.78 
64.01 


0.469 
0.467 
0.465 
0.462 
0.457 


2.0 
1.8 
1.6 
1.4 


0.4 
0.4 
0.4 
0.4 


1.85 
1.85 
1.85 
1.85 


10.0 
10.0 
10.0 
10.0 


i.24 x ur 7 
1.72 x ur 7 
2.46 x Hr 7 
3.67 x ur 7 


0.0757 
0.0923 
0.114 
0.144 


1.97 x 10~ 4 
2.48 x 10~ 4 
3.23 x 10~ 4 
4.44 x 10~ 4 


75.49 
75.21 
72.92 
69.14 


234.20 
198.88 
167.40 
138.92 


0.492 
0.491 
0.490 
0.489 



The 1.4 Mq trajectory in Fig. 1 1 with R v oc r 1 / 3 , for L%\ = 3.4. 



Table 3 

Modified Evolutionary Wind Models: = 3.4 and Yf ~ 0.46 



Mass (Mq) 


L™(10 51 ergs" 1 ) 


R w (km) 


M(M Q s" 1 ) 


2(10 48 ergs" 1 ) 


/VchflO^ergs- 1 ) 


t p (ms) 




1.8 


15.7 


9.0 


1.01 x 10~ 5 


6.18 


0.092 


5.63 


146.44 


1.6 


15.7 


9.0 


1.48 x 10~ 5 


7.79 


0.122 


5.61 


122.98 


1.4 


15.7 


9.0 


2.29 x 10~ 5 


10.2 


0.174 


5.39 


102.15 


1.4 


15.7 


8.5 


1.94 x 10~ 5 


9.25 


0.157 


4.97 


108.48 


1.4 


15.7 


8.0 


1.62 x 10~ 5 


8.35 


0.141 


4.54 


115.86 
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FIG. 1. — Matter velocity (v) in cm s 1 as measured in the Schwarzschild frame as a function of radius (R) in km for a 1.4 Mq (gravitational) protoneutron 
star with L,-, e = 8, 7, 6, 5, 4, 3, 2, and lxlO 51 erg s _1 . For L„ e = 8 X 10 51 erg s'\ we set (e Ve ) = 11 MeV, (e De ) = 14 MeV, and (e Vfl ) = 23 MeV. For each 

subsequent luminosity, the average neutrino energy for each species was decreased according to (e„} oc L l J 4 . The luminosities were set in the ratios LvjL Ve = 1.3 
and L,y e /L, y(J = 1 .4. The dots mark the critical point for each wind profile, where the matter velocity is equal to the local speed of sound. The neutrinosphere radius 
is held fixed at 10 km. For all profiles, for r < 60 km, the flow is nearly homologous with v oc L v r. 
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FIG. 2. — Logio of the mass density (p) in g cirr 3 versus radius (R) in km for the same range of neutrino luminosities as in Fig. 1 and for the same protoneutron 
star characteristics. Dots mark the critical point. The inset shows log 10 p versus R for the region close to the protoneutron star. Note the steep density gradient, which 
drops precipitously over as much as five orders of magnitude in just over a kilometer. As the neutrino luminosity decreases, p(R v ) increases in order to maintain our 
inner integral boundary condition on the v e neutrino optical depth (eq. 18). 
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FIG. 3. — Matter temperature (T) in MeV versus radius in km for the same profiles as in Fig. 1. Dots mark the critical point. Note that the important regime of 
possible r-processing lies between 0.5 MeV and approximately 0.08 MeV, a region extending out to 700 km for the highest luminosities and to less than 200 km for 
the lowest luminosities shown. Comparing the temperature in this range of radii to those in Fig. 1 , it's clear that the assumption of constant outflow velocity (Meyer 
and Brown 1997) is a poor approximation during nucleosynthesis in neutrino-driven winds. The inset shows the structure of the temperature as a function of radius 
very close to the neutrinosphere. The bump in T is caused by the onset of heating (compare with Fig. 4). 
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FIG. 4. — Energy deposition rate (q) in units of 10 21 erg g 1 s 1 as a function of radius (in km) for the wind models in Figs. 1-3. q profiles for v e neutrino 
luminosities from 8 X 10 51 erg s -1 to 1 X 10 51 erg s -1 are depicted. The total heating rates (g; eq. 14) for each of these models are given in Table 1. 
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FIG . 5 . — Entropy (s, per baryon per ks ) versus radius in km for the same protoneutron star wind models as in Figs. 1 -4. Note that the entropy asymptotes quickly; 
comparing this figure to Fig. 3, we can see that in most cases * increases less than 10 units for T < 0.5 MeV. 
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FIG. 6. — Contributions of specific neutrino processes to the energy deposition rate for = 8 X 10 51 erg s -1 , R v = 10 km, and M = 1.4 Mq. The solid line shows 
the total heating rate due to all processes, the short-dashed line is the net contribution from the charged-current reactions (q cc ) (eqs. 20 and 21), the long dashed line 
is neutrino-electron/positron scattering (q ve \ eq. 27), the short dot-dashed line is the net energy deposition rate due to i/,P; <-» e + e~ (q U D\ eqs. 29 and 31), and the 
long dot-dashed line is for neutrino-nucleon scattering (q V N \ eq. 28). Note the fairly rapid decrease in q cc at r ~ 35 km is due to the recombination of free neutrons 
and protons into a particles. As the neutrino luminosity decreases in these models, this transition moves inward in radius so that, for the lowest luminosities, the 
charged current processes end abruptly at ~ 25 km. For models with larger M/R v q ue and q V v become more important relative to q cc . However, even for M = 2.0 
Mq and R v = 10 km, q cc dominates at the peak in <j t ot- 
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FIG. 7. — The electron fraction (Y e , solid line), proton fraction (X p , short dashed line), neutron fraction (X„, long dashed line), and alpha fraction (X a , dot-dashed 
line) for the highest luminosity wind model (Lp e = 8 X 10 51 erg s~') in Figs. 1-5, corresponding to the heating profile in Fig. 6. Note the transition at r ~ 35 
km where a particles form and effectively shut off the charged-current heating rate. The asymptotic Y e is set very close to the neutron star (r ~ 15 km) and only 
undergoes small subsequent change as a result of the a-effect (Fuller and Meyer 1995; McLaughlin, Fuller, and Wilson 1996)). In this case, Y e changes by ~1% 
beyond r = 20 km. The asymptotic electron fraction, Y", is well-approximated by eq. (16) in the text. 
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FIG. 8. — Dynamical timescale (t p ) defined in eq. (32) versus radius for the eight models shown in Figs. 1-5. On each line, the first point marks the radius at 
which T(r) = 0.5 MeV. The second point marks where T(r) = 0. 1 MeV. Note that for high luminosities, r p increases by more than a factor of three over the range 
of radii and temperatures relevant for r-process nucleosynthesis. For the lowest luminosities dr p /dr changes sign and r p actually decreases by ~30% between the 
two points on the L De = 1 X 10 51 erg s -1 curve. 
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FIG. 9. — Tracks for constant protoneutron star mass in the plane of asymptotic entropy (s a , per ks per baryon) and dynamical timescale (r p , eq. 32) in seconds. 
Ri/ is held constant at 10 kilometers. Each track covers a range in luminosity from La e = 8.0 X 10 51 erg s _1 to Lc e = 0.7 X 10 5 ' erg s _1 . Small dots are the points 
for which a model with a given luminosity was calculated. All masses are gravitational masses. In these models, we take {e Ve ) = 11 MeV, (ea e ) = 14 MeV, and 

(e Vfl ) = 23 MeV for L„ e = 8.0 X 10 51 erg s~'. For each subsequent luminosity, the average energies were decreased according to (e v ) oc L l J A . In Table 1, we 
summarize the global properties at several luminosities along the 1.4 Mq trajectory in this figure. 
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FIG. 10. — A schematic showing total neutrino luminosity in 10 51 erg s -1 (thin solid line), average v e neutrino energy inMeV (long dashed line), and neutrinosphere 
radius (R„) in km. Two possible evolutions for R v (t) are shown. Linear contraction such that R v (t = 0.4s) = 20.3 km and R v (t = 1 s) = 10.0 km, which we label as 
'R v (t) oc 1 -at', is shown as a short dashed line. Contraction with R v (t) oc f~'/ 3 is shown as a dot-dashed line. The thick solid line denotes 10 km, the final R v for 
the protoneutron star. L'° l is proportional to r 09 . We set (s^^)/{eo e ) = 1.6, {ev e )/{e Ve ) = 1.3, LjjjL v ^ = 1.4, and LjjjL Ve = 1.3. 
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FIG. 1 1 . — Evolutionary wind models in the plane of asymptotic entropy (.«„ per ks per baryon) and dynamical timescale (r p , eq. 32) in seconds. Tracks with 
Ri/(t) oc (1 — at) (see Fig. 10) for M = 1.4, 1.6, 1.8, and 2.0 Mg are shown as solid lines with small dots. The solid line without dots is the track with R v (t) oc f~'/ 3 
for 1.4 Mq. All tracks start at small s a and r p with R v ~ 20.3 km and L|' = 8.0. The tracks end at high s a and long r p with L S p = 0.4. Taking L v oc f~° 9 , this 
range spans the first ~ 12 seconds of protoneutron star wind evolution. The large open circles markLj' = 3.4 for each mass trajectory. Dashed lines connecting 
all mass tracks with large dots are lines of constant log 10 [AM c j] (see eq. 34, Fig. 13, and text for details). The dashed line connecting all mass tracks on the right, 
marked with -6.0, delineates at what point in the evolution only 10"* M of material is yet to be ejected. In all of the more than 350 models presented here the 
asymptotic electron fraction is always in the range 0.45 < ¥" < 0.495. Generally, Y" increases as the luminosity and average neutrino energies decrease. Finally, 
the dash-dotted lines in the upper left corner indicate, for constant values of Y", in what range of s and r p production of the third r-process abundance peak is 
possible in the calculations of Meyer & Brown (1997). See §6 for discussion and details. 
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FIG. 12. — The log of the integrated total mass ejected M e j (eq. 34) in units of Mq as a function of time in seconds, assuming L„(f) tct oc r - 9 for the tracks of 
constant mass shown in Fig. 11. Evolutionary models with oc (1 — at) (see Fig. 10) are shown as solid lines with small dots, which correspond to the dots in 
Fig. 1 1, indicating luminosity (time). M e ; for the track in Fig. 1 1 with R u {t) oc ?~'/ 3 and M =1.4 Mq is shown here as a long dashed line without dots. The large 
dot on each of the solid curves at ( ~ 1 second shows the point at which the model has contracted to R v = 10 km. This corresponds to the point on Fig. 1 1 where a 
given track of constant mass with R v (t) ct(\-at) takes a sharp turn at r p ~ 0.006 seconds. 
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FIG. 13. — Logio[AM c j(=M^j"—M e j (())], in units of Mq, versus time in seconds for the models in Fig. 12. The long dashed line is for the model with R v (t) act '/ 3 . 
The four solid lines with dots correspond to those in Figs. 11 and 12 for the 1.4, 1.6, 1.8, and 2.0 Mq evolutionary models, which we label with 'R v {t) oc 1 -at'. 
Lines of constant log 10 [AM e j] are shown in Fig. 1 1 as dashed lines connecting big dots on the evolutionary models with Rv{t) oc (1 -at). 



